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We measure the mass of the top quark in lepton+jets final states using the full sample of pp collision data 
collected by the DO experiment in Run II of the Fermilab Tevatron Collider at yT = 1.96 TeV, corresponding to 
9.7 fb _1 of integrated luminosity. We use a matrix element technique that calculates the probabilities for each 
event to result from tt production or background. The overall jet energy scale is constrained in situ by the mass 
of the W boson. We measure m, = 174.98 ±0.76 GeV. This constitutes the most precise single measurement of 
the top-quark mass. 
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I. INTRODUCTION 
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Since the discovery of the top quark IH 12, the determi¬ 
nation of its properties has been one of the main goals of 
the Fermilab Tevatron Collider, and more recently also of the 
CERN Large Hadron Collider (LHC). The measurement of 
the top-quark mass m t , a fundamental parameter of the stan- 
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dard model (SM), has received particular attention, since m t , 
the mass of the W boson, Mw, and the mass of the Higgs boson 
are related through radiative corrections that provide an inter¬ 
nal consistency check of the SM 0. Furthermore, m, dom¬ 
inantly affects the stability of the SM Higgs potential, which 
has cosmological implications mm 

The world-average m, = 173.34 ±0.76 GeV has reached 
a precision of about 0.5% 00, which has been limited for 
some time by systematic uncertainties. The main systematic 
uncertainties arise from the calibration of the jet energy scale 
(JES) and the modeling of tt events in Monte Carlo (MC) sim¬ 
ulations. 

In this manuscript, we present a measurement of m, based 
on the full set of pp data at y/s = 1.96 TeV recorded by the 
DO detector in Run II of the Fermilab Tevatron Collider, cor¬ 
responding to an integrated luminosity of 9.7 tb 1 . The anal¬ 
ysis focuses on tt events identified in /'-l-jets final states with 
t representing either an electron or a muon. This includes 
contributions from leptonic x —t £V/V r decays. The top and 
antitop quark are assumed to always decay into a W boson 
and a b quark ED- Thus, at leading order (LO) in perturba¬ 
tive quantum chromodynamics (QCD), one of the W bosons 
in the W + W~bb final state decays into a charged lepton and 
a neutrino, while the other decays into a pair of quarks, and 
all four quarks ( qq'bb ) evolve into jets. Such £±jets events 
are characterized by an isolated electron or muon with large 
transverse momentum pj, a large imbalance in transverse mo¬ 
mentum caused by the undetected neutrino, and four high -p? 
jets. We exploit this distinct signature in the event selection to 
discriminate tt events from background. 

To extract the value of m t , we utilize a matrix element 
(ME) technique that determines the probability of observing 
an event under both the tt signal and background hypotheses 
described by the relevant ME ifTTIl . The overall JES is cal¬ 
ibrated in situ through a constant factor fcjt.s by taking into 
account all the kinematic information in a given event, while 
constraining the reconstructed invariant mass of the W —t qq' 
decay products to Mw = 80.4 GeV fl2l . 

We provide a detailed description of the measurement that 
was presented previously in a Letter (El, and is an update 
and supersedes a previous DO measurement of m, = 174.94 ± 
1.14(stat±JES) ±0.96(syst) GeV ll4l using 3.6 ftV 1 of in¬ 
tegrated luminosity. In the present measurement, we use not 
only a larger data sample to improve the statistical precision, 
but also refine the estimation of systematic uncertainties, re¬ 
sulting in a reduction of the uncertainty by 50%. This im¬ 
provement is achieved through an updated detector calibra¬ 
tion, in particular through improvements to b quark JES cor¬ 
rections na, and through the use of recent improvements in 
modeling the tt signal. In the following, we point out the most 
important improvements of this analysis with respect to the 
previous measurement. The analysis is performed “blinded”, 
i.e., with no reference to the resulting m, , until the final ap¬ 
proval of the analysis methodology. 

As in all direct measurements of m t , this analysis relies on 
simulated tt MC events, which are used for the absolute cal¬ 


ibration of the response of the ME technique. Therefore, the 
extracted m t corresponds to the generated mf en parameter im¬ 
plemented in the LO MC event generator used to simulate tt 
events, which in our case is ALPGEN m. Although mf en is 
not well defined at LO, it is expected to correspond within 
« 1 GeV to the m, defined in the pole-mass scheme 03 . 
There is currently no universal MC event generator available 
for simulating tt events at next-to-leading order (NLO) that 
includes the effects of parton showering and hadronization. 
Programs such as MC@NLO m or POWHEG ED , commonly 
referred to as NLO generators, simulate only the qq —> tt , 
gg —t tt, and qg —> tt processes at NLO, whereas the decays of 
the top quark are simulated at LO through the narrow-width 
approximation. 


This manuscript is arranged as follows. A brief description 
of the DO detector is given in Section[EI] which is followed in 
Section |III]by a summary of simulations used to model signal 
and background events. The selection of tt candidate events is 
reviewed and the sample composition presented in Section[IV] 
Section[V]discusses the ME technique used to extract in ,, with 
particular emphasis placed on the improvements since the pre¬ 
vious publication 03. This is followed in Section [Vl]by the 
description of the measurement of the signal fraction f in data 
using the ME technique, including its calibration. The mea¬ 
surement of m, and the calibration of the response of the ME 
technique in m, is discussed in Section |VII| After describing 
the evaluation of systematic uncertainties in Section [VIII| we 
compare our measurement with previous results in Section|IX] 
and conclude and give the final result with statistical and sys¬ 
tematic uncertainties in Section [X] Appendix |A| discusses the 
issue of initial and final state radiation, and Appendix |B] sum¬ 
marizes the parametrization of the detector response through 
transfer functions. 


II. T HE DO DETECTOR 


The DO detector (20ll2Tl l contains a magnetic central track¬ 
ing system, calorimetry, and a muon system. The central 
tracking system comprises a silicon microstrip tracker (SMT) 
and a central fiber tracker (CFT), both located within a 2 T 
superconducting solenoidal magnet. The SMT ll22l has « 
800000 individual strips, with typical pitch of 50 — 80 gm, 
and a design optimized for track and vertex finding within 
1 77 1 < 2.5 l23l . The system has a six-barrel longitudinal struc¬ 
ture, each with a set of four layers arranged axially around 
the beam pipe, and interspersed with radial disks. In 2006, a 
fifth layer, referred to as Layer 0, was installed close to the 
beam pipe |24|| . The CFT has eight thin coaxial barrels, each 
supporting two doublets of overlapping scintillating fibers of 
0.835 mm diameter, one doublet being parallel to the colli¬ 
sion axis, and the other alternating by ±3° relative to the 
axis. Light signals are transferred via clear fibers to solid- 
state visible-light photon counters (VLPCs) that have ss 80% 
quantum efficiency. 

Central and forward preshower detectors, located just out- 
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side of the superconducting coil (in front of the calorimetry), 
are constructed of several layers of extruded triangular scintil¬ 
lator strips that are read out using wavelength-shifting fibers 
and VLPCs. These detectors provide initial sampling of elec¬ 
tromagnetic showers, and thereby help distinguish between 
incident electrons and jets. The next layer of detection in¬ 
volves three liquid-argon/uranium calorimeters: a central sec¬ 
tion (CC) covering up to |rj | ss 1.1, and two end calorimeters 
(EC) that extend coverage to |rj| ss 4.2, housed in separate 
cryostats. The electromagnetic (EM) section of the calorime¬ 
ter is segmented into four layers, with transverse segmenta¬ 
tion of the cells in pseuodorapidity and azimuth of A/] x A© = 
0.1 x 0.1, except for the third layer, where the segmentation is 
0.05 x 0.05. The hadronic portion of the calorimeter is located 
after the EM sections and consists of fine hadron-sampling 
layers, followed by more coarse hadronic layers. In addition, 
scintillators between the CC and EC cryostats provide sam¬ 
pling of developing showers for 1.1 < 1 77 1 < 1.4. 

A muon system ll25l is located beyond the calorimetry, and 
consists of a layer of tracking detectors and scintillation trig¬ 
ger counters before 1.8 T iron toroid magnets, followed by 
two similar layers after the toroids. Tracking for 1 77 1 < 1 re¬ 
lies on 10 cm wide drift tubes, while 1 cm mini-drift tubes are 
used for 1 < I'll <2. 

Luminosity is measured using plastic scintillator arrays lo¬ 
cated in front of the EC cryostats, covering 2.7 < |T 7 1 < 4.4. 
The trigger and data acquisition systems are designed to ac¬ 
commodate the high instantaneous luminosities of the Teva- 
tron ||20l [26l . Based on coarse information from tracking, 
calorimetry, and muon systems, the output of the first level 
of the trigger is used to limit the rate for accepted events to « 
2 kHz. At the next trigger stage, with more refined informa¬ 
tion, the rate is reduced further to « 1 kHz. These first two 
levels of triggering rely mainly on hardware and firmware. 
The third and final level of the trigger, with access to all of 
the event information, uses software algorithms and a com¬ 
puting farm, and reduces the output rate to ss 100 Hz, which 
is recorded. 


III. THE DATA SAMPLE AND SIMULATIONS 


The sample of pp collision data considered in this anal¬ 
ysis was collected between April 2002 and September 2011 
in Run II of the Fermilab Tevatron Collider. The sample is 
split into four data-taking epochs: “Run Ha”, “Run Ilbl”, 
“Run IIb2”, and “Run IIb3” corresponding to 1.1 tb 1 , 
1.2 ftT 1 , 3.0 lb 1 , and 4.4 tb 1 of integrated luminosity, re¬ 
spectively, after implementation of data-quality requirements. 
All MC simulations are split accordingly into epochs to model 
variations in the response of the detector, such as the inser¬ 
tion of Layer 0 into the SMT detector between Run Ha and 
Run Ilbl El, and the dependence of the track reconstruction 
efficiency on instantaneous luminosity l27l . 

Simulated MC events for tt signal and the dominant W+jets 
background, which are used to verify the agreement be¬ 


tween data and simulations and to calibrate the measurement 
(described in Sections | VI A| and VII A[ >, are generated us¬ 
ing the MC generator ALPGEN m to simulate the hard- 
scattering process with up to five additional partons, and 
PYTHIA ll28l l29l to simulate hadronization and shower evo¬ 
lution. The MLM matching scheme l30l is employed to avoid 
overlaps between components of the event belonging to the 
hard process, implemented through a ME, and parton evolu¬ 
tion (showering) into jets. The W +jets background samples 
are divided into two categories: (i) VV' -I hf which represents 
W +cc and W +bb production, possibly in association with u, 
d, s quarks, or gluons, and (ii) W +lf which represents all other 
flavor configurations. The cross section for the W 4 jets contri¬ 
bution provided by ALPGEN is multiplied by a factor of 1.3 to 
match the NLO calculations of the MCFM MC generator ED- 
In addition, the heavy-flavor contributions of category (i) are 
multiplied by a factor of 1.47, using the same calculations. 

We also consider other small background contributions on 
the order of few percent to improve the agreement between 
data and simulations. The Z+jets background is simulated in 
the same way as W +jets with ALPGEN, taking into account 
contributions from Z+hf production, with rates normalized to 
NLO calculations ED. PYTHIA is used to simulate diboson 
(WW, WZ, or ZZ) production, while COMPHEP+PYTHIA ll32l 
is used for single top quark processes that are simulated 
for m, = 172.5 GeV. The hard-scattering process for all 
signal and background MC samples is simulated using the 
CTEQ6L1 set of parton distribution functions (PDF), except 
for COMPHEP, which uses the NLO CTEQ6M set 13^1 . More 
details on the MC simulations and their normalization can be 
found in Ref. 


The simulation of parton showers with PYTHIA uses tune A 
in combination with the CTEQ6L1 PDF set and with Aqcd = 
0.26 GeV, also referred to as “DO tune A”. The detector re¬ 
sponse is fully simulated through GEANT3 ESI- To simulate 
the effects from additional pp interactions (pileup), events 
taken with minimal trigger requirements, selected from ran¬ 
dom pp crossings at the same instantaneous luminosity as the 
data, are overlaid on the fully simulated MC events. These 
events are then reconstructed using the same algorithms as 
used on data. 


Events from multijet (MJ) production can pass our selection 
criteria, which happens typically when a jet mimics an elec¬ 
tron or when a muon arising from semileptonic decay of a b or 
c quark appears isolated. The kinematic distributions from the 
MJ background are modeled using data events for which lep¬ 
ton isolation requirements are inverted. The absolute contri¬ 
bution of this background to each of the e+jets and /iljets fi¬ 
nal states is estimated using the “matrix method,” described in 
Ref. EH- This method uses the number of events with leptons 
originating from the hard interaction, ’ , and from MJ pro¬ 
duction, JVjm , both with less restrictive lepton-identification 
criteria, and relates them to the contributions to the sample of 
events that fulfil standard lepton identification criteria. This 
is calculated using the relation N = £"' W N[^ s + e MJ N^ se , 
where e ,uW and e MJ represent the efficiency of events that pass 







6 


TABLE 1: A summary of kinematic event selections. 


One charged lepton 

pj > 20 GeV 

Ini < 1-1 (e) 

|tj| <2.0 (p) 

Four jets 

pj > 20 GeV 

Ini < 2-5 

Jet with highest pj 

pj > 40 GeV 

Ini <2.5 

Imbalance in transverse momentum 

jj T > 20 GeV 



the loosened lepton identification criteria to also pass the stan¬ 
dard identification criteria. The parameters " and e MJ are 
measured, respectively, in control regions dominated by lep¬ 
tons originating from the hard interaction and by MJ events. 


IV. EVENT SELECTION AND SAMPLE COMPOSITION 


The trigger selects f f jets events by requiring either at least 
one lepton (electron or muon), or at least one lepton and at 
least one jet, resulting in an efficiency of 95% for tt candidate 
events containing an electron and 80% for events containing a 
muon, before any offline selections are applied. 

Selected tt candidate events are required to have at least one 
primary pp collision vertex (PV) with three or more tracks re¬ 
constructed within 60 cm of the center of the detector in the 
coordinate along the beam axis. In addition, there must be 
exactly four jets with pj > 20 GeV within 1 77 1 < 2.5, with 
the leading jet further satisfying pj > 40 GeV. Jets are re¬ 
constructed using an iterative cone algorithm (36l with a cone 
parameter of R cone = 0.5. Jet energies are corrected to the 
particle level using calibrations derived from exclusive y+jet, 
Z+jet, and dijet events fi~5l . These calibrations account for 
differences in detector response to jets originating from a 
gluon, a b quark, or u.d.s. and c quarks. Furthermore, we re¬ 
quire the presence of one isolated electron ED or muon m 
with transverse momentum pj > 20 GeV and /] < 1.1 or 
1 77 1 <2, respectively. The electrons and muons are required to 
originate from within 1 cm of the PV in the coordinate along 
the beam axis and to be separated in A R = \J (A 77) 2 + (A (j)) 2 
from the closest jet by A R > 0.5, where At] is the difference 
in pseudorapidity and A (j) in azimuth (in rad) between the 
lepton and the jet. At least one neutrino is expected in the 
£+jets final state, and an imbalance in transverse momentum 
of ]jj > 20 GeV is required. The /)-/■ variable is defined as the 
opposite of the vector sum of the transverse energies in each 
calorimeter cell, corrected for the energy carried by identified 
muons and after applying the jet energy scale calibration. The 
kinematic selections are summarized in Table |T] 

Events containing an additional isolated charged lepton, 
i.e., an electron with pj > 10 GeV and |r]| < 2.5 or a muon 
with pj > 10 GeV and |tj <2, are rejected. In addition, 
any fi +jets events with an invariant mass /V/^, of the isolated 
muon and a second loosely isolated muon with pj > 15 GeV 
within 70 < <110 GeV are rejected to suppress Z+jets 

events. 

To reject fi +jets events where the momentum of the muon 


has been mismeasured, we require \jj < 250 GeV and Mj < 
250 GeV, where Mj is the transverse mass of the W bo¬ 
son l38l , To further reduce contributions from such events, 
which is found at the % level, we use additional requirements 
on the significance of the curvature of the muon track, S K . It is 
defined as the curvature K = l/pr divided by its uncertainty, 
where Q is the electric charge of the muon. We reject events 
with muons within one or both of the two regions where they 
are likely to be mismeasured, S K < A(j)(fi,j/j) x 25.5 — 70 or 
S K < A(j)(fi,rfj) x 4.4 + 8 . 8 , where 7 -) is the difference 
in azimuth between the muon momentum and the direction of 
1 It- 

To reduce the contribution from MJ production caused by 
the misidentification of a quark or gluon jet as a lepton, which 
changes the calculated jfj, we require A <j)(l,rfj) > 0.5. Addi¬ 
tional reduction of the MJ background is achieved by requir¬ 
ing A(j)(e,^j) > 2.2 — ifj x 0.045 GeV -1 in the e+jets final 
state and A(j)(p,rfj) > 2.1 — rfj x 0.035 GeV~* in the p+ jets 
final state. 


Since two Zxquark jets are expected in the final state, at least 
one jet per event is required to be tagged as originating from 
a b quark (fc-tagged) through the use of a multivariate algo¬ 
rithm |[39l . The tagging efficiency in this analysis is ss 65% 
for /xquark jets, while the mistag rate for gluons and for u.d.s 
quark jets is se 5%. 


We select 1502 events in the e+jets final state in data, 
while 1286 are selected in the /i+jets final state. To verify 
that the simulations agree with the data following all selec¬ 
tions, we use the NNLO calculation, including next-to-next- 
to-leading-logarithmic (NNLL) corrections |40l . to normal¬ 
ize the tt signal. The contribution from W +jets events is 
adjusted to normalize the absolute yield from our signal and 
background models to data before applying b jet identification 
criteria. The relative yields of contributions from W+(bb. nip) 
and W + (cc,«lp) production and contributions from W+n lp 
production are further adjusted using data, as described in 
Ref. 1341 . The number of observed events is compared to 
expectations in Table [2] considering statistical uncertainties 
only. 


The distribution in the invariant mass m t f of the tt sys¬ 
tem is shown in Fig. [26] (a) and (b) in Section |VIIC to¬ 
gether with predictions from our signal and background mod¬ 
els. The distribution is generated with a mass of the top quark 
of mf cn = 172.5 GeV. The mass m t j is constructed using 
j)t as an initial estimate for pj of the neutrino, assuming 
Mw = 80.4 GeV, and solving the resulting quadratic equa¬ 
tion for pi. When no solution exists, p v r is scaled to obtain a 
solution. The resulting neutrino momentum is then combined 
with the momenta of the four jets and the charged lepton to 
obtain m t j. The invariant mass, M\v, of the two jets that match 
one of the W bosons is presented in Fig. [27] (a) and (b) in 

is unique 


Section VIIC The reconstruction of Mw in Fig. 27 


for events with two /?-tagged jets, as the two other jets are as¬ 
signed to the W boson. If only one jet is ^-tagged, all possible 
assignments of jets to partons in the LO approximation of tt 
decays are tried, and the one yielding the smallest value of 
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TABLE 2: Comparison of the number of observed data events after 
all selections with the expectation from simulations, assuming a t j = 
7.24 pb~' (40) and using mf en = 172.5 GeV. The category “other 
backgrounds” encompasses Z+jets, WW, WZ, ZZ, and single top 
quark production, as well as dileptonic it decays. The production of 
a W boson in association with at least one b or c quark is denoted as 
“W + hf’, while “W + If’ stands for all other flavor configurations. 
The uncertainties are purely statistical. 


Contribution 

e+jets 

p+jets 

tt 

918.1 

± 

3.6 

824.9 

± 

3.5 

Other backgrounds 

97.8 

± 

0.5 

79.2 

± 

0.9 

Vk + hf 

126.0 

± 

2.1 

162.2 

± 

2.8 

W + If 

77.9 

± 

2.1 

101.0 

± 

2.9 

Multijet 

144.4 

± 

24.2 

48.2 

± 

16.1 

Expected 

1364.1 

± 

24.7 

1215.5 

± 

17.0 

Observed 


1502 



1286 



|nif —m, I is retained. If the /?-tagged jet is matched to the top 
quark decaying into three jets, the two untagged jets matched 
to the same top quark form the W boson. If all three jets that 
are matched to this top quark are not /;-tagged, the two low¬ 
est pr jets are matched to the W boson. A similar match¬ 
ing algorithm is used if there are more than two ^-tagged 
jets. The uncertainties in the predictions include the dominant 
sources of systematic uncertainty that change the shape of the 
kinematic distributions: JES, jet energy resolution, jet iden¬ 
tification, flavor dependence of the detector response, b tag¬ 
ging, lepton identification, lepton momentum scale, trigger, as 
well as hadronization and higher-order effects for simulated tt 
events. Good agreement is observed between data and simula¬ 
tions. We verify that the agreement improves using the values 
for the signal fraction /, m t , and kjgs measured in this analy¬ 
sis, as shown in Section|VIIC| 


V. ANALYSIS METHOD 

The extraction of m t is based on the full kinematic and topo¬ 
logical information in the event, and is performed through 
a likelihood technique using probability densities (PDs) for 
each event, described by the ME of the processes contributing 
to the observed events. 


A. Matrix element technique 

The ME technique was first suggested by Kondo ED, and 
was experimentally applied for the first time at the Tevatron 
in the measurement of m t CD, and later in the determination 
of the helicity of the W boson ll42ll and in the first evidence 
obtained for production of single top quarks mm. 

Assuming two non-interfering contributing processes, tt 


and W +jets production, the per-event PD is: 

■Pevt = A(x)\fP (x; m t , £jes) 

+ (1 -/)fl>kg(*;&JEs)] J (1) 

where the observed parameters /, m t , and Zcjes, are to be de¬ 
termined from data, x HU represents a vector of the mea¬ 
sured four-momenta of the four jets and the charged lepton, 
but not ifj due to its limited experimental resolution, and .4 (x) 
accounts for acceptance and efficiencies. The function Psig 
describes the PD for tt production, while Pbkg describes the 
PD for W +jets production. W +jets events, i.e., the sum of 
W + hf and W + If categories in Table [2j represent the domi¬ 
nant background contribution of 14% in the e+jets and 20% 
in the ti+jets final states, derived with the normalization pro¬ 
cedure described in Section UV] The next-to-dominant back¬ 
ground contribution is from MJ events. It is accounted for 
in Pevt through F\, k a , since W +jets and MJ production have 
similar PDs in the selected kinematic region. The similar¬ 
ity of PDs for W 4 -jets and MJ production can be understood 
in that the main difference between the two processes arises 
from the presence of the W boson in the W +jets production, 
while an energetic jet which is misidentified as an isolated 
lepton is present in MJ production. At the same time, both 
processes are characterized by the presence of four additional 
jets, which, together with the W boson for W +jets production 
or the jet misidentified as an isolated lepton for MJ production, 
set a similar momentum-transfer scale of the hard interaction 
process. The effect of including MJ production through Tj,kg 
is accounted for by the calibration procedure outlined in Sec¬ 
tions [VI] and IVIII The combined contribution from all other 
backgrounds amounts to about 6 % in both final states (cf. Ta¬ 
ble [ 2 ), and is therefore not included in the calculation of P ev t- 


B. The signal probability 

The probability density for tt production to yield a given set 
of partonic final-state four-momenta y l45l in the hard scatter¬ 
ing of two massless quarks with four-momenta q\ and <72 is 
proportional to the differential cross section dcr of the corre¬ 
sponding process, given by 

. (2n)* \JLq^ t t\ 2 

d<Ttf-wf(?;m,) = ’ 1 ft' J ' dT> 6 , ( 2 ) 

2g M v<7i‘72 

where ./Mqq^u denotes the ME for the qq—yti^t b(£v)b(qq') 
process, g^y is the tensor of the Minkowski metric, and d<t >6 
is an infinitesimal element of six-body phase space. 

To obtain the differential cross section in pp collisions, the 
expression in Eq. (|2} is convoluted with the corresponding PD 
for all possible flavor combinations of the colliding quark and 
antiquark, yielding 

d (J p p-> t t(y;m t ) =Y< d qi d q2f{q\)f{q2) 

quark _ 
flavors #1 > <72 

X dOqq-Ht(y m , m t ). (3) 
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The longitudinal momentum PDFs, f(qi, z ), are taken from 
CTEQ6L1 ll33l . while the dependence of /(<;/ J X ) and f{qi. y ) 
on transverse momentum is taken from the PD in the PYTHIA 
simulation. Only quark-antiquark annihilation at LO is taken 
into account in -//qq-,tt, and the differential cross section 
d Opp^n therefore does not represent the full differential cross 
section for tt production in pp collisions. Effects from gluon- 
gluon and quark-gluon induced tt production, accounting for 
about 15% of the cross section, as well as from higher-order 
corrections from real emissions, are accounted for in the cali¬ 


bration procedure described in Sections VIA and VIIA 


In general, the set x of measured four-momenta is not iden¬ 
tical to the set of corresponding partonic variables y because 
of finite detector resolution and the evolution of quarks into 
jets. As will be discussed in Section V D| this is taken into 
account through the transfer function (TF) W (x,y,k]Es), rep¬ 
resenting the probability density for the measured set x to have 
arisen from the partonic set y. Thus, the probability to observe 
a given reconstructed tt event characterized by x is obtained 
through a convolution with the TF in the calculation of the 
differential cross section 


do p p^ t T(x;m t ,kjEs) = dy da pP ^ t j(y\m t ) 

Jy 

x W(x,y\k]Es) ■ (4) 


The PD to observe a tt event with a set of kinematic quan¬ 
tities x observed in the detector is given by 


1. Calculation of 


The LO matrix element for the qq —> tt process 061071 is 
given by 

\.^qq^a\ 2 = [6n g ^ FF- (2-^ 2 sin 2 (0 gf )) , (8) 

where a s is the strong coupling, /3 = \p t \/E t represents the ve¬ 
locity of the t (or t) quark in the qq rest frame, and d qt denotes 
the angle between the incoming parton and the outgoing top 
quark in the qq rest frame. The form factors FF are identical 
to those given in Eqs. (24) and (25) of Ref. 08]. 


2. Calculation of G pp ^ t t,obs 


The calculation of the total cross section for tt produc¬ 
tion observed in the detector defined in Eq. 0 is challeng¬ 
ing, because a double integral has to be performed over all 
possible partonic and measured configurations for each set of 
(m t ,k]Es)- We therefore apply a MC integration technique, 
and factorize the partonic and measured parts of the calcula¬ 
tion to obtain: 


o P p-+tT,obs{"h) = dydcj p p^ t f(y;?n t ,k 3 Es)A(x)W(x,y-,k 3 Es) 

Jy 


X (4 {ttl; , kjEs))i 


(9) 


where 


P s ig(x;/n,,kj ES ) 


4 6 (a, ui tl k jes) 

tfpp—>7?,obs 5 &JES ) 


(5) 


where the total cross section for tt production observed in the 
detector, defined as 


Gpp - r It. 



J X 


Gpp^tt 


(x; m,,k]Es), 


( 6 ) 


Gpp—ttt, tot 



dy d Gpp-ytt 


(f, m t ,kjEs) 


( 10 ) 


is the total partonic cross section for tt production in pp col¬ 
lisions as shown in Fig. [T] and (A(m t ,kms))x I s the mean ac¬ 
ceptance averaged over all possible configurations of x for a 
given set of (m t ,kjEs) parameters. 


3. Calculation of the mean acceptance 


ensures a proper normalization of A(x)P s ; g . Using Eqs. (j2jl, 

0- and Q, the PD to observe a tt event with four-momenta x The calculation of the mean acceptance (. A(m t , &jes))x P ro_ 

in the detector can be written explicitly as ceeds in two steps: 


Psig(x;m t ,kjEs) = 


x 


x 

x 


1 


Gpp—±tt, obs ip^t j &JES ) 

E f dqidq 2 f{qi)f{q 2 ) 

quark _ 
flavors quqi 


(27T) 4 | 

'igpvq^q'i 


2 

-dd> 6 


W(x,y k S Es) ■ 


(7) 


The parametrization of \^f q q->tt\ 2 an< l ^pp-ttt, obs are discussed 


below, and the calculation of P sl „ is given in Section V F 


(i) First, (A(m,,kjEs = l))jis determined using MC events 
processed with the full simulation of the detector, ac¬ 
cording to 


(A(m,,kms = !))•? 


L Nacc 

i =T ^MC^exp 

L ^Vgen 

i= i WMC 


(ii) 


In this expression, V gen denotes the total number of gen¬ 
erated events, and A' acc represents the total number of 
events from within detector acceptance, while wmc is 
the MC weight from ALPGEN, and vr exp the weight that 
accounts for experimental factors such as the trigger and 
identification efficiencies. 
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m t [GeV] 



m, [GeV] m t [GeV] 



m, [GeV] m, [GeV] 


FIG. 1: The total partonic cross section for tt production in pp 
collisions through quark-antiquark annihilation at LO, as defined in 
Eq. {TO}. 


(ii) In the second step, the dependence of the mean accep¬ 
tance on A;jes is derived using simulated parton level 
events, where the detector response is taken into ac¬ 
count through the TF W(x,y; Tjes) for the correspond¬ 
ing kjEs value. 

The mean acceptance is shown for each data-taking epoch in 
Fig.[2]for e+jets and in Fig. [3] for fi 4-jets final states. 


C. The background probability 


The differential partonic cross section for 1]^ is calculated 
according to 

-Pbk g (r;£jEs) = —-or—7 

®pp— >VK+jets,obs V^JES ) 

X E / dc tUzdq 2 .zf(qi,z)f{q2,z) 

quark ^ 
flavors ^l, z ’^2,z 

^ (IK) | tS^qq— >VK+jets | 

2gpv<fiq2 6 

x W(x.y\ kjEs) , (12) 

which is similar to P s ; g in Eq. {7}, with the exception that the 
initial-state partons are all assumed to have p j = 0. The ME 
jets is calculated at LO using the ME for IT+4 jets 
production as implemented in the VECBOS |49] program. 
Since VECBOS treats the final state quarks as massless, the 
same ME is used to model the production of W bosons in as¬ 
sociation with u,d,s,c, and T-quark jets. By definition, Tkkg 
is independent of m,. The dependence of l]^„ on A'n :s is taken 
into account, as before, through the TF. 


FIG. 2: Mean acceptance (A(«7 ,,Ajes = l))r with respect to the 
total partonic cross section for pp —>■ tt production, as a function of 
m t for different £jes values in the e+jets final state, for (a) Run Ha, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. The acceptance tends to 
decrease for later epochs due to a decreased reconstruction efficiency, 
which is mainly caused by the higher instantaneous luminosity in the 
second half of Run II and detector aging. 



m, [GeV] m t [GeV] 



m, [GeV] m t [GeV] 


FIG. 3: Same as Fig. [2] but for the /i+jets final state for (a) Run Ha, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. 


D. Modeling of detector response 


Here, we describe the modeling of the detector response 
through the TF W(x,y; Ties)- We use a parameterized TF, as 
a full simulation of the detector would not be computationally 
feasible. 
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In constructing the TF, we assume that the functions for in¬ 
dividual final-state particles are not correlated. We therefore 
factorize the TF into contributions from the four jets and the 
charged lepton used to calculate P s { g . The resolution of the 
jfj is worse than that of the p-y of jets and leptons, and is 
therefore not used in defining event probabilities. We assume 
that the directions of the e, jJ ., and jets in (rj, 0 ) space are 
well-measured. The effect of finite detector resolutions in 77 
and <j) is studied for jets and found to be negligible. The TF 
for jet and charged lepton directions are therefore represented 
by 8 functions, S 2 (r), (j)) = S(ri y — T] x )S(</) y — <j) x ), thereby re¬ 
ducing the number of integration variables by 5x2 = 10 . 
The finite resolutions on the energies of jets and charged lep¬ 
tons are explicitly taken into account, as described in Sec¬ 
tions rVD~nrVD2l and |VD3l 

There is an inherent ambiguity in assigning measured jets 
to partons from tt decay. Consequently, all 24 jet-parton as¬ 
signments that are possible in LO are considered. The in¬ 
clusion of b tagging information provides improved identifi¬ 
cation of the correct jet-parton assignment through weights. 
The weight w, for the i-th jet-parton assignment is defined 
by a product of individual weights wj for each jet j. For 
b tagged jets, wj is equal to the per-jet tagging efficiency 
fc'tag ( CX k ; Ej. rjj), where CXy labels the configurations where 
the jet is assigned to either (i) a b quark, (ii) a c quark, or 
(iii) a u,d,s quark or gluon. For untagged jets, the wj fac¬ 
tors are equal to 1 — $a g (a*; E- r . rj 7 ). Hence, the weights \v t 
are close to unity for those configurations that assign /?-tagged 
jets to b quarks and untagged jets to u.d.s quarks or gluons, 
and vv, <C 1 otherwise. Since the contributions to W +jets 
processes are parameterized by ^^^w+jets without regard to 
heavy-flavor content, the weights vv,- for each of the 24 permu¬ 
tations in the background probability are all equal. The sum 
of weights is always normalized to unity: 1 w i = 1 - 

We define the TF as 


W{x,y; kms) = W e (E x ,Ey)Si(r],^) 

24 f 4 j 

x L w ‘ j n (13) 


(=1 


J =1 


where E is the energy and W; et is the TF appropriate for the en¬ 


ergy resolution of jets, as discussed in Section V D 1 and Wy 
is either W e , the energy resolution for electrons, as discussed 
in Section |V D 2| or VV),, the resolution in the transverse mo¬ 


mentum of muons, as discussed in Section [V D 3| 


1. Modeling of jet momentum 


The TF for jets represents the probability that the jet energy 
E x measured in the detector corresponds to a parent quark of 
energy E y . It is parameterized in terms of a double Gaussian 
function whose means and widths depend on E y . For Icjes = 1, 


it is given by 


W ]et (E x ,E y -k m s = 1) 


1 

V2tz(P2 + P3P5) 


X 



+ exp 


[(E x -E y )- Pl } 2 \ 

2p\ ) 

[{E x -E y )-p A ] 2 \ 

2 Pi J 


(14) 


where the pi are linear functions of the energy of the parent 
quark: 


Pi =a i + bi E y . 

For kjEs 7 ^ 1, the jet TF changes to 


Wjet (E x ,E y ; kjEs) = 


«fc(4’ E r ;1 

&JES 


(15) 


(16) 


where the kjES factor in the denominator preserves the nor¬ 
malization JdE x Wj et (E x ,E y ;kjEs) = T 

The parameters a, and b,- in Eq. are determined from 
fully simulated tt events for each data-taking epoch, after 
applying the jet energy corrections and resolution smearing, 
which are required to match the detector response in data. 
To avoid potential bias from the kinematic selections, all the 


requirements described in Section IV are applied at detec¬ 
tor level, except for having exactly four jets, at least one of 
which is ^-tagged, and except the p-y threshold for the jets, 
which is lowered to 15 GeV. The parameters a,- and bj are 
extracted from MC generated tt events with m t values rang¬ 
ing between 150 and 190 GeV, for each data-taking epoch. 
We use an unbinned likelihood fit that maximizes the prod¬ 
uct of the Wjet terms for all jets in all events as a function 
of at and b L . Separate sets of a,- and b, parameters are de¬ 
rived for three categories of jets: (i) from light quarks (m, d, 
s ) and charm, (ii) from b quarks with a muon of a relatively 
low pr within the jet cone (soft muon tag), and (iii) for all 
other b quarks. This is done for four regions of 77 : 1 77 1 < 0.4, 
0.4 < 1 77 1 < 0.8, 0.8 < 1 77 1 < 1.6, and 1.6 < 1 77 1 < 2.5 to take 
into account the 77 dependence of the detector response. These 
same four regions are used to derive the correction for the 
flavor-dependence of the detector response to jets fI51 . When 
deriving the TF, the matching of jets to partons is performed 
using a simple cone-matching algorithm, where a jet is consid¬ 
ered uniquely matched to a parton if A/v'fjet. closest parton) < 
f^cone/2 = 0.25. The cone-matching parameter is chosen to 
maximize the matching efficiency, while minimizing the mis¬ 
matching rate of quarks to jets from pile up. The probability 
of a detector-level jet not to be matched to a particle-level jet 
is 2% per jet. Figure[4]illustrates the TF for jets from category 
(i) as a function of E x for different values of E y for Run IIb2, 
which is representative of all four data-taking epochs. The TF 
for jets from category (ii) is shown in Fig. [5] and for cate¬ 
gory (iii) in Fig. [ 6 ] The values of the a,- and parameters for 
each of the four data-taking epochs are summarized in Ap¬ 
pendix [b] 
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FIG. 4: The TF M / j et (£ t ,£,,;/:jp_s) for light-quark jets at kjps = 1 in 
Run IIb2, as a function of the measured jet energy E x , for different 
parton energies E y , as indicated in the legend. The TF are shown in 
four pseudorapidity regions: (a) |rj| < 0.4, (b) 0.4 < 1 77 1 < 0.8, (c) 
0.8 < 1 77 1 < 1.6, and (d) 1.6 < |tj| < 2.5. Energies E x > 150 GeV are 
also considered, but not shown. 



FIG. 6: Same as Fig. [4] but for />quark jets without muons inside 
the jet cone. The TF are shown in four pseudorapidity regions: (a) 
|rj| < 0.4, (b) 0.4 < |tj| < 0.8, (c) 0.8 < 1 77 1 < 1.6, and (d) 1.6 < 
M<2.5. 


Fig ,[7]we compare the distributions in the invariant mass of the 
dijet system that matches the two quarks from the W —7 qq' de¬ 
cay, and the invariant mass of those two quarks, smeared with 
the jet TF. We compare the distributions in the invariant mass 
of the trijet system matched to one of the top quarks and the 
corresponding quarks smeared with the jet TF in Fig. [8] Both 
tests indicate that the jet TF well describe the jet resolutions. 
Any residual disagreements are taken into account through the 


calibration procedures, described in Sections VIA and VIIA 


2. Modeling of electron momentum 


The electron energy resolution is parameterized by the TF 


W e {E x ,Ey) 


-exp 

\f2no 


1 ( E x -E> 

2\ a 


21 


(17) 


with 


FIG. 5: Same as Fig. [4] but for fo-quark jets with at least one muon 
inside the jet cone. The TF are shown in four pseudorapidity regions: 
(a) 1 77 1 < 0.4, (b) 0.4 < 1 77 1 < 0.8, (c) 0.8 < 1 77 1 < 1.6, and (d) 1.6 < 
M<2.5. 


To verify that the TF describes the detector response to jets, 
we use fully simulated ffMC events with all jet energy correc¬ 
tions and smearing needed to match resolutions in data. They 
must pass all selections summarized in Section IV and each 
of the four jets needs to be uniquely matched to a quark. In 


E' 


= E y + 0.324 GeV, 


<7 = 


S = 


y/(0.028 • E'y) 2 + (S ■ E'y) 2 + (0.4 GeV) 2 , 


0.164 GeV 2 0.122 GeV 






- exp 


C 


sin 9 e 


C = 1.3519- 


2.0956 GeV 6.9858 GeV z 


E 'y 


E? 


(18) 
(19) 

-C, (20) 

( 21 ) 


where E x is the reconstructed electron energy, E y is the elec¬ 
tron energy at parton level, and 0 e is the polar angle of the 
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FIG. 7: Invariant mass of the dijet system matched to the W —> qq' 
decay for (a) Run Ha, (b) Run Ilbl, (cj Run IIb2, and (d) Run IIb3. 
Fully simulated jet-parton matched MC events passing all selection 
requirements are compared to partons with energies smeared using 
the jet TF, and then applying the same standard event selection crite¬ 
ria. 
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FIG. 8: Same as Fig. [7] but for the invariant mass of the trijet sys¬ 
tem matched to the hadronic decay of the top quark for (a) Run Ha, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. 


electron relative to the proton beam direction. These param¬ 
eters are obtained from the modeling of electron energy re¬ 
sponse and resolution in Ref. lf50l for Run Ha, and are used 
for all four data-taking epochs because they changed only 
marginally between the epochs. 


TABLE 3: Parameters of the muon TF for muons with at least three 
associated hits in the SMT and for all other tracks. These values are 
the same for all the four data-taking epochs. 


Parameter 

> 3 hits 
in the SMT 

< 3 hits 
in the SMT 

so (GeV- 1 ) 

2.082 x 10" 3 

3.620 x 10- 3 

■«1 

1.125 x 10- 2 

1.388 x 10- 2 

c 0 (GeV- 1 ) 

7.668 x 10 -3 

2.070 x 10- 2 

Cl 

7.851 x IQ- 2 

7.042 x IQ- 2 


3. Modeling of muon momentum 


The resolution of the central tracker, which provides the 
measurement of the muon pj, is described through the un¬ 
certainty on the curvature K = \/pt of the muon track. The 
muon TF is parameterized as 


w 'n (K x ,Ky) 


1 

= J — ex p 

\J7mg 


1 (K*-Ky 

2\ a 


2 ' 


( 22 ) 


where K x is the reconstructed muon curvature, and K y is the 
corresponding quantity at parton level. 

The resolutions 

f s for |t]| < 1.4 | 

| y // s 2 + {c-(|?7|-1 .4)} 2 for |tj| > 1.4 J (23) 

are obtained from muon tracks in simulated events, where the 
s and c parameters are linear functions of K r : 

s = i 0 +si K y , (24) 

c = co + ci Ky. (25) 

The values of these coefficients, which are the same for all 
four data-taking epochs are given in Table [3] for muon tracks 
with at least three associated hits in the SMT, and for all other 
tracks. This simplified parameterization of the momentum 
resolution is valid at high transverse momenta (pj > 20 GeV) 
where the limitations in coordinate resolution dominate over 
the effects of multiple scattering. 


E. Extraction of the parameters of interest 


In this Section, we discuss how the parameters /, m t , 
and kjEs are extracted from a sample of N selected events 
that are characterized by their measured quantities X,v = 
{xi,X 2 , ...,xn}- A likelihood function, «Sf(X^;m r ,kjEs,/), is 
constructed at each grid point from the product of the individ¬ 
ual P e vt values according to 


N 

m t X JES,/) = n^evtK; m ti ^JES 5 f ) • (26) 

i= 1 
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To construct this likelihood from P e vt, /Tig is calculated 
on a grid in (m t ,k jes) with spacings of (1 GeV,0.01). For 
simulated tt events, it covers 25 points in m, ranging be¬ 


tween mf — 12 GeV and mf + 12 GeV, and 21 points 
in kjES between —0.1 and kffff, +0.1. The ranges are 
chosen to ensure that they contain the numerically relevant 
part of the likelihood for all pseudo-experiments (cf. Sec¬ 
tions |VI A| and |VII A). For data, the grid is extended to 
[153 GeV, 192 GeV] in m, and to [0.85,1.15] in k ms . This 
choice is made because the values of m, and Ties i n data are 
not known a priori , as the analysis is performed blinded, fol¬ 
lowing the procedure to be described in Section |VIIB For 
simulated background events, the same grid is chosen as for 
data. Since Fj,kg does not depend on m ,, it is calculated on a 
grid in £jes only, with the same spacing and ranges as P s ; g . 


1. Extraction of the signal fraction 


The likelihood is determined by the value of the signal frac¬ 
tion f(m t ,kjEs) that maximizes f£ for every assumed pair of 
(m r ,kjEs) values from the grid on which 2z? is calculated. As 
will be detailed in Sections [VI] and |VII| the measurement of 
the signal fraction is performed prior to the measurement of 
m, and kjES- To extract / fit independently, the parameters m, 
and kjEs are integrated over a uniform prior 

,.fit _ /dWfdkjEs /.Sf(Xjv; m t ,k JES ,f) 
f dw,dkjEs ^ (X^; m ,, kj E s, /) 


by summing over all points (m r ,kjEs) on the grid that reflect 


the numerically relevant part of the likelihood, cf. Section V E 


2. Extraction of the top-quark mass and kj£g 


To obtain the best unbiased estimate of mf 1 , 
the two-dimensional likelihood 2z?(Xjv; m t ,£jes) = 
2z^(Xjv; m tl k]ES:f{" l t 1 k]Es)) is projected onto the m, 
axis 

-S?(Xa r, m,) = y dkjEs 2 z?(Xa r, m,,kjEs) (28) 

using Simpson’s rule ED- The best unbiased estimate for mf 1 
is then given by 


fit f dm, m,) 

ft, = - 

Jdm, _S?(Xjv; m t ) 


and its uncertainty by 


cr fit = 

mt 


f dm, (in, — mf 1 ) 1 m t ) ^ 1 

f dm, 2Sf(Xjv; m,) 


(29) 


(30) 


Thus, by construction, the statistical uncertainty on m, in¬ 
cludes the impact of the statistical uncertainty on £jes- To 
obtain the best unbiased estimates for kf^ s and <7*^, the same 
procedure is applied as to estimate mf 1 and by replacing 


/tjes <->• m, in Eqs. (28][30i. 


F. Calculation of the signal probability 


In the calculation of P sl „ in Eq. |7|, a total of 24 variables are 
associated with the integration over the momenta of the two 
colliding quarks dtj\dqn and over the phase space of the six- 
body final state, d<t>6. After taking into account the 8 2 (ri,(j)) 
functions for the directions of the four jets and the charged 
lepton according to the TF of Eq. the dimensionality 
of the integral is reduced by 2 x 5 = 10. Imposing energy- 
momentum conservation provides 4 additional constraints, re¬ 
sulting in a 10-dimensional integral that must be evaluated. 
The transverse momentum of the neutrinos is inferred from 
momentum conservation, resulting in up to eight solutions for 
the neutrino kinematics, which are averaged. The measure¬ 
ment of /)-/ is not used because of its limited experimental 
resolution. 

The integration in Eq. 0 is performed numerically using 
MC integration 1152} . To reduce the computational demands 
of the integration, importance sampling {53} is used. To fully 
exploit importance sampling, we perform a Jacobian transfor¬ 
mation of the nominal 10 integration variables to optimized 
variables where prior information is either known or can be 
obtained easily. This prior information is then used in the im¬ 
portance sampling. Thus, the integration is performed over the 
masses of the top and antitop quarks, which are assumed to be 
equal, the masses of the W + and bosons, the energy (cur¬ 
vature) of the electron (muon), E q /(E q +Eg,) for the quarks 
from the W —t qq decay in the LO approximation, and the 
transverse momenta < 7 i jX ,< 7 i,y,< 72 ,x,g 2 ,y of the colliding quarks. 
The constraint M w = 80.4 GeV for the in-situ JES calibration 
is imposed by integrating over Mw according to a prior given 
by a Breit-Wigner distribution. More details about the moti¬ 
vation for the choice of integration variables can be found in 
Ref. |54). 

The computation of P s j g in the previous version of this anal¬ 
ysis M, which is based on about 1/4 of the full Run II inte¬ 
grated luminosity, and uses only one set of MC simulations, 
required two million CPU-hours {55). To make our ME tech¬ 
nique applicable to the full Run II data sample of 9.7 tb 1 , 
and to perform a calibration of the method with dedicated MC 
simulations for each of the four data-taking epochs, we imple¬ 
ment two new approaches to reduce computational demand. 
These approaches are described in detail in Ref. l54} . and are 
briefly summarized below. 

(i) We utilize low-discrepancy sequences for the numerical 
MC integration of Fsig according to Eq. ([7]), instead of 
conventional pseudo-random numbers in order to bene¬ 
fit from the key property of low-discrepancy sequences 
of sampling the integration phase space in a maximally 
uniform way. This leads to a convergence rate of the es¬ 
timated integral inversely proportional to /V s;imp | e (56), 
which is superior to y/AI samp i e expected for pseudo¬ 
random numbers {53), where A+, mp i e is the number of 
samplings of the integrand. We utilize the implementa¬ 
tion of Bratley and Fox {57) of the Sobol sequence (58) . 
which is among the best performing low-discrepancy 
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sequences. 

An accurate, and computationally efficient estimate of 
the numerical precision of the integral is a key ingre¬ 
dient of numerical MC integration. Since the standard 
error estimator based on the variance of the integrand is 
too pessimistic for low-discrepancy sequences, because 
of their superior convergence rate, we implement a ded¬ 
icated estimator of the numerical uncertainty according 
to Ref. l59l . which represents a novel treatment in the 
context of the ME technique. 

The utilization of low-discrepancy sequences is of par¬ 
ticular interest, as it is a procedure universally applica¬ 
ble to MC integration and independent of the computing 
hardware. Their implementation in our ME technique 
provides a reduction of about one order of magnitude in 
computation time, compared to pseudo-random num¬ 
bers, for the same numerical accuracy. 


(ii) We factorize the £jes parameter from the ME compu¬ 
tation, a novelty in the context the measurements of m t 


with the ME technique. As described in Section V E 


Aig has to be calculated for 25 x 21 = 525 grid points 
in (m tl kjEs) for simulated tt events. Instead of recalcu¬ 
lating P s ; g for each point in we calculate the 

integral 


E f dqidq 2 f(qi)f(q2) ^ d<f> 6 (31) 

luark _ J _ A£jUV#l#2 

iavors QhQl 

in Eq. (jT|i only once for a given in ,, and use the result 
to calculate P s ; g for all kjES,<=i,2,...,2i through multipli¬ 
cation by W(x,y;kjEs,!')- This yields another order of 
magnitude reduction in computation time. 


In total, we obtain a reduction of computation time from about 
2 h per event, averaged over the sample of simulated tt events 
for m, = 172.5 GeV, to about 80 s per event, i.e., by about two 
orders of magnitude. Both improvements are verified Il54l to 
provide a performance of the ME technique consistent with 
that in Ref. Ifl4l . They prove essential for a dramatic improve¬ 
ment in the precision of this analysis by allowing to increase 
the number of simulated MC events used for the calibration of 
the ME technique (cf. Sections | VI A| and | VII A|), and for the 
evaluation of systematic uncertainties (cf. Section|VIII|i. 


VI. MEASUREMENT OF THE SIGNAL FRACTION 


The ME technique calculates P e vt from first principles with 
a LO ME and a parameterized description of the detector re¬ 
sponse through W(x,y;kjEs)- The calculation of the parame¬ 
ters of interest, in this case /, has to be therefore calibrated 
through simulations. The calibration of the response of the 
ME technique in / and the subsequent extraction of / from 
data are the subject of this Section. 
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FIG. 9: The response of the ME technique in /, derived from the 
average of / fit over 1000 PEs at each / gen point, for (a) Run Ha, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3 in the e+jets final state. 
The response is fitted with a linear function shown as a solid black 
line, whose slope and offset values are given at the bottom of the 
plots. The broken red line represents an ideal response. 


A. Calibration of the method 


Fully simulated tt events, as described in Section III at 
mf en = 172.5 GeV and kfy'^ = 1, expected from the standard 
JES calibration procedure |jT31 , are used for calibrating the re¬ 
sponse of the ME technique in /. Similarly, fully simulated 
W +jets events at /fjj™ = 1 are included to model the dominant 
background contribution from W (jets production. We refine 
the calibration procedure relative to that of Ref. Ifl4l by explic¬ 
itly accounting for the next-to-dominant contribution from MJ 
production by using a sample of events in data obtained as de¬ 
scribed in Section III Since the detector response, and there¬ 
fore its parametrization, changes between data-taking epochs 
and final states, this calibration is obtained for each epoch and 
final state. 


At each of the five generated / gen points, ranging between 
0.5 and 0.9 in steps of 0.1, 1000 pseudo-experiments (PEs) are 
constructed, each with the same number of events as observed 
in data. Since the matrix method provides a normalization 
of the background contribution from MJ production, each PE 
contains the same number of randomly drawn MJ events as 
determined through the matrix method, as shown in Table [2] 
The remaining number of background events is then drawn 
from the sample of simulated W -Fjets events, according to a 
binomial distribution in 1 —/ gen . 

The PE ensembles at each of the five / gen points are used 
to obtain a linear calibration for the response of the ME tech¬ 
nique in /. This is achieved by parametrizing the signal frac- 
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FIG. 10: Same as Fig. [9] but in the /i+jets final state, for (a) Run Ila, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. 


tion / fit , extracted as for data in SectionlVEland averaged over 


the PEs for each given / gen point, versus / gen . The resulting 
linear parameterizations of the response of the ME technique 
in / are shown for each data-taking epoch in Fig. [9] for the 
e+jets, and in Fig. [TO] for the /i+jets final states. Typical 
slope parameters are close to 0.75, while typical offsets are 


few percent at / gen ss 0.70, as expected from Section IV 


B. Results 


The extracted / fit in data is calibrated using the results from 
the previous Section for each final state and epoch according 
to 

/•fit _ 

f=- - L , (32) 

s f 

where oj and Sf are the offset and slope parameters, respec¬ 
tively, as given in Figs. |9]and|T0] The uncertainties on oj and 
sf are propagated to / assuming they are Gaussian distributed. 
The resulting / are given in Table [4j split by final state and 
by epoch. Averaging / over the entire Run II, weighting the 
contributions from data-taking epochs by their respective inte¬ 
grated luminosities, yields / = 63% in the e+jets and / = 70% 
in the /i+jets final states, with an absolute uncertainty of 1% 
from the finite size of the data sample and the calibration. 

As a cross-check we translate the measured / values into 
cross sections for tt production, which are also reported in Ta¬ 
ble]?] Combining them in a similar fashion to /, we find a t j = 
7.8 ±0.1 (stat) pb in the e+jets and <7„- = 7.6 ±0.1 (stat) pb 
in the /i+jets final states. The cross sections determined with 
the ME technique are in good agreement with the prediction 


TABLE 4: Signal fractions / measured in data after calibration, split 
by final state, and for each data-taking epoch as well as combined 
over the entire Run II. Also shown are the resulting measured cross 
sections for tt production. The combined uncertainty from limited 
statistical accuracy in data and from the finite number of simulated 
events used for calibration amounts to 1% for /. 


Epoch 

Final state 

Signal fraction 

<7„- (pb) 

Run Ila 

e+jets 

0.72 

8.9 

/i+jets 

0.65 

7.8 

Run Ilbl 

e+jets 

0.77 

7.6 

/i+jets 

0.66 

6.8 

Run IIb2 

e+jets 

0.68 

7.8 

/i+jets 

0.66 

7.5 

Run IIb3 

e+jets 

0.56 

7.6 

/i+jets 

0.75 

8.0 

Run II 

e+jets 

0.63 

7.8 

/i+jets 

0.70 

7.6 


of <7 „- = 7.24^q 2 y pb (scale+PDF) from a NNLO calculation 
at m t = 172.5 GeV, including NNLL corrections [|40l , and us¬ 
ing the MSTW2008NNLO set of PDFs gO). This is also in 
agreement with o t j = 7.78+q^J pb measured by the DO Col¬ 
laboration in 5.3 fb 1 of integrated luminosity ED- 


VII. MEASUREMENT OF THE TOP-QUARK MASS 


As described in Section [yT] the response of the ME tech¬ 
nique in m t and kjEs has to be calibrated prior to the mea¬ 
surement of those parameters. In this Section, we describe the 
calibration procedure, report the measured m t and kjES values, 
demonstrate that simulations with our measured values of /, 
m t , and ±n;s describe the data well, and draw comparisons to 
other measurements of m t . 


A. Calibration of the method 


The calibration of the response of the ME technique in m t 
and Ges proceeds similarly to the calibration in f described in 
Section [VT] It uses seven simulated samples of tt events, five 
at mf n = 165,170, 172.5,175,180 GeV for = 1, and two 
at = 0.95,1.05 for mf en = 172.5 GeV. These are com¬ 
plemented by three simulated samples of W +jets events at 
^■jes = 0-95,1, and 1.05. The calibration procedure is refined 
relative to Ref. m by accounting explicitly for the next-to- 
dominant contribution from MJ production. Together, the tt, 
W +jets, and MJ samples are used to obtain a linear calibra¬ 
tion for the response of the ME technique in m t and kjEs ■ The 
calibration is applied separately for each data-taking epoch 
and final state to account for varying detector response and 
sample composition. 

At each of the generated points, 1000 PEs are 
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FIG. 11: The response of the ME technique in m t , obtained using 
simulations from the average of mf x in 1000 PEs at each mf en point, 
as a function of mf e , for (a) Run Ha, (b) Run Ilbl, (c) Run IIb2, and 
(d) Run IIb3 in the e+jets final state. The response is fitted with a 
linear function shown as a solid black line, whose slope and offset 
values are summarized at the bottom of the plots. The red broken 
line represents an ideal response. 


constructed, each containing the same number of events as 
observed in data, by randomly drawing simulated signal and 
background events according to the measured / for each data- 
taking epoch and final state, cf. Table [4] For each PE, the 
fraction of signal events is changed randomly according to a 
binomial distribution around the / value measured in data. 
Since the background conUibution from MJ production is 
modeled using data, it is not changed for different generated 
(mf en , ) points. 

The response of the ME technique in m, and kj^s is ob¬ 
tained by extracting mf l and kSL for each pseudo-experiment 
as i n dat a, i.e., according to the procedure described in Sec¬ 
tion VE The averages (mf 1 ) and (kf^ s ) are determined from 
Gaussian fits to distributions in mf 1 and kf^ s , respectively. A 
linear parametrization of the response in m, is obtained by fit¬ 
ting (ml' 1 ) as a function of mf en . Those parameterizations are 
shown for each epoch in Fig. 11 for the e+jets and in Fig. [12] 
for the ll +jets final states, together with the offset and slope 
of the linear fits. The response in m, is close to the ideal case 
of no offset and unit slope. At mf en = 172.5 GeV, which is 
close to the current world average of m t = 173.34 GeV 0, 
the difference between the obtained and ideal response in m t 
is ~ 0.5 GeV, i.e., at the level of « 0.25% - a very good 
performance given that a LO ME with parameterized detec¬ 
tor response is used in the ME calculation. We verify that 
the ideal response of the ME technique is retained through 
ensemble studies using parton-level tt events generated with 
PYTHfA using a LO ME. In a similar fashion, we obtain a lin- 


FIG. 12: Same as Fig. El but in the /r+jets final state, for 
(a) Run Ha, (b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. 


ear parametrization of the response of the ME technique in 
kjEs, as shown in Fig. 13 for the e+jets and in Fig. [14] for the 
ji +jets final states. As for the case of m t , the response is close 
to ideal. At = 1 closest to the expectation from the stan¬ 
dard JES correction ED, the difference between the obtained 
and ideal response is at the level of 1%. 


To check the validity of the ct,„ ; and <7/, |rs uncertainties ex¬ 
tracted from the ME estimates described in Section |V E| we 
calculate the pulls in m t (n mi ) and in kjgs (% JES X according to 


Tim, - 


nat™ — 


mf 1 — (mf 1 ) 


ufA 

k jes 


rrfit 

^m t 

~(k] 


fit \ 
JES/ 


7 fi* 

*-'JES 


(33) 

(34) 


For an ideal response, the pull should follow a Gaussian dis¬ 
tribution of unit width centered at zero. In our case, both 7Z mi 
and ^ JES are Gaussian distributed, and centered within uncer¬ 
tainties at zero. For each of the simulated (m,,kjEs) points, 
we determine the pull widths from the a parameters of Gaus¬ 
sian fits to n tltl and 7R. jes . The average width of n mt is then 
obtained from a constant fit of the distribution of pull width 
, as shown for each epoch in Fig. 15 for e+jets 


versus m\ 


gen 


and in Fig. 16 for /i+jets final states. The average width of 
7t,t JES is obtained in the same way and is given for each epoch 
in Fig. [FT] for e+jets and in Fig. 18 for /I + jets final states. 
Typically, pull widths are ss 20% higher than unity, indicat¬ 
ing that the uncertainties cr,^ and <s estimated with the ME 
technique are somewhat underestimated. They are therefore 
corrected by inflating by the value of the pull width found for 
the respective epoch and final state. 
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FIG. 13: The response of the ME technique in &jes> obtained using 
simulations from the average of in 1000 PEs at each point, 
as a function of for (a) Run Ha, (b) Run llbl, fc) Run IIb2, 
and (d) Run IIb3 in the e+jets final state. The response is fitted with 
a linear function shown as a solid black line, whose slope and offset 
values are summarized at the bottom of the plots. The red broken line 
represents an ideal response. The poor compatibility of data with the 
linear fit displayed in panel (a) likely to occur if many fits, in our 
case, 16, are performed. 


FIG. 15: The width of the distribution of the pull in m t , obtained 
using simulations at each m, point, as a function of mf en , for 
(a) Run Ha, (b) Run llbl, (c) Run IIb2. and (d) Run IIb3 in the e+jets 
final state. The error bars are given by the uncertainty on the width 
of the distribution in K nh , which is not corrected for the effect that the 
same event can appear in many PEs, and therefore somewhat under¬ 
estimated. The pull is fitted with a constant shown as a solid black 
line, whose value is given at the bottom of the plots. The red broken 
line represents an ideal pull width of unity. 
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FIG. 17: The width of the distribution of the pull in £jes> th¬ 
rived using simulations at each point, as a function of k||g, for 
(a) Run Ha, (b) Run llbl, (c) Run IIb2. and (d) Run IIb3 in the e+jets 
final state. The error bars are given by the uncertainty on the width of 
the distribution in 7t ^, which is not corrected for the effect that the 
same event can appear in many PEs, and therefore somewhat under¬ 
estimated. The pull is fitted with a constant shown as a solid black 
line, whose value is given at the bottom of the plots. The broken line 
represents the case of ideal pull width of unity. 
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FIG. 18: Same as Fig. |T5] but in the p+jets final state, as derived 
for (a) Run lla, (b) Run llbl, (c) Run IIb2, and (d) Run IIb3. 


B. Results 


Here, we present the results of our application of the ME 
technique to the full Run II data sample corresponding to 
9.7 fb 1 of integrated luminosity. The analysis is performed 
blinded in m t by adding a constant unknown offset to the 
extracted m t value. The offset is drawn according to a uni¬ 
form distribution in the range [—2 GeV,2 GeV]. This in¬ 
terval roughly corresponds to [—3<7 wa, 3(7wa], where g\va = 
0.76 GeV denotes the total uncertainty on the current world 
average in m t . 


The method to extract m t and kjES is described in Sec¬ 
tion |VE The linear calibrations from Section [YII A are ap¬ 
plied to data through the likelihoods J 2 ?(Xjv; m t ), defined in 
Eq. (28 1 , and 2z?(X,/v; £jes) by transforming the m, and kj^s 
parameters according to 


, {m t - 172.5 GeV) — o m 
m, —y m t = -1- 172.5 GeV,(35) 


kms —■> ^jes — 


(^JES ~ 1) “OAtjes 


s kn 


+ 1 j 


(36) 


where the offset and slope parameters are given in Figs. m 
[14] The direct calibration of the likelihoods 22?(X,v; m t ) and 
^jes) takes into account the effect of slopes ± 1 when 
estimating the statistical uncertainties a mt and (T/ || s . Subse¬ 
quently, the slope-corrected uncertainties (7„ J; are calibrated 
by multiplying them with the average widths of the distribu¬ 


tions in n mt , summarized in Fig. 15 for e+jets and in Fig. 16 


for /i+jets final states. Similarly, the uncertainties <Jk ms are 


calibrated with the average widths in 7R JES , given in Figs. 17 
andfl8l 


The m, measurements are summarized in Fig. 19 along 
with their statistical uncertainties, which include the statistical 
contribution from /cjhs- The observed statistical uncertainty 
is compared with expectations from pseudo-experiments for 

for e+jets and in Fig. 


gen 


= 172.5 GeV in Fig. 


20 


21 


for 


ll +jets final states. The observed uncertainties are consistent 
with expectations from simulation. We obtain the combined 
m, for the full Run II data set of 9.7 ftG 1 by calculating the 
uncertainty-weighted mean, which is the best linear unbiased 
estimator in case of Gaussian uncertainties, according to 


L 5 rn t,i 

1=1 jj 7 

(m t ) = 


1 




<7 = 


1 ®m t ,i, 


(37) 


where i runs over all final states and epochs. The combined 
result after unblinding is m, = 174.98 ± 0.58 (stat+IES) GeV. 
To verify that the values measured in each epoch and final 
state are consistent within their statistical (stat+JES) uncer¬ 
tainties, we calculate the value of the '/} distribution for 7 
degrees of freedom, which yields a probability of 62%. 


The measured kjES values are shown in Fig. 22 together 
with their statistical uncertainties. We also compare the ob¬ 
served values of 07 JES with expectations from simulation, and 
find consistency. We obtain a combined value of &jes = 
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FIG. 19: Summary of the measured m, values with their statistical 
uncertainties for each data-taking epoch and split according to final 
state. The red-shaded band indicates the statistical uncertainty on 
the combined m, value in Run II data corresponding to 9.7 fb -1 of 
integrated luminosity. 
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FIG. 21: Same as Fig. [20] in fl+ jets final states in (a) Run Ha, 
(b) Run Ilbl, (c) Run IIb2, and (d) Run IIb3. 
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FIG. 22: Same as Fig. |T9] but for the £jes parameter. 


FIG. 20: Distribution of expected statistical uncertainties G nh in 
e+jets final states in (a) Run Ha, (b) Run Ilbl, (c) Run IIb2, and 
(d) Run IIb3, obtained in 1000 PEs at mf en = 172.5 GeV and kj^s = 
1, after calibrating for non-unit width of the distribution in Tt„ h . The 
observed statistical uncertainties in data are indicated by arrows. 


1.025 ± 0.005 (stat+Wf) by applying Eq. 37 after replacing 
kjES tn t . The individual measured kjES are consistent with a 
X 2 probability of 54%. The total JES uncertainty for tt events 
after all selections is between 2.0% and 2.1%, depending on 
the epoch, consistent with our measured £jes at the level of 


1.3 standard deviations (SD). Due to the busy environment of 
tt events containing four jets, a higher fraction of the event’s 
energy is expected in jet cones leading to a higher £jes ac¬ 
cording to Eq. (fl6|), which is not the case in y+jet and dijet 
samples where the JES is derived. 


Figure [23] presents the two-dimensional likelihood density 
in (/WfjkjEs), which is obtained by multiplying the calibrated 
likelihood densities for each data-taking epoch and each final 
state. We compare the m, and £jes values obtained from the 
minimum of the double-Gaussian fit to the two-dimensional 
likelihood density in Fig. 23 and find them to match at the 
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FIG. 23: Two-dimensional likelihood (X#; m/,^jEs)/-^max for 
data. Fitted contours of equal probability are overlaid as solid lines. 
The maximum is marked with a cross. Note that the bin boundaries 
do not necessarily correspond to the grid points on which «Sf is cal¬ 
culated. 
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FIG. 24: (a) Distribution of the expected uncertainty on m t , as 

obtained using PEs where the tt contribution is taken at m, = 
172.5 GeV and £jes = 1, after calibration. The observed statisti¬ 
cal uncertainty is indicated by an arrow, (b) Same for uncertainty on 
fc JES- 




level of 10 4 . In Fig. |24lwe compare the observed statistical 


uncertainties on m t and £jes with expectations from PE studies 
using tt samples for mf en = 172.5 GeV and /cj^ = 1. The 
observed uncertainty is consistent with the expectation. 

The values of in, and £jes measured in <?±jets and ti+jets 
final states using 9.7 fb 1 of integrated luminosity are 


mf** 

, c+jets 

'tjES 

g+jets 

m, 

,g+j e ts 

k jes 


175.55 ±0.81 (stat+JES) GeV, 
1.026 ±0.006 (stat+m f ), 

174.36 ± 0.84 (stat+JES) GeV, 
1.025 ±0.007 (stat+m,). 


The corresponding two-dimensional likelihood densities in 
(m t ,kms) and the comparison between the expected statistical 
uncertainty for each final state and the observation are shown 
in Fig. [25] The results in e±jets and /i±jets channels are 
consistent at the level of 1 SD considering just statistical un- 


FIG. 25: (a) Two-dimensional likelihood .if(XAr; m f ,£jEs)/-^max 

for data in the e+jets final state. Fitted contours of equal proba¬ 
bility are overlaid as solid lines. The maximum is marked with a 
cross, (b) Distribution in expected uncertainty on m, in e+jets fi¬ 
nal states, obtained using PEs where the tt contribution is taken at 
mf en = 172.5 GeV and £jes = 1, after calibration. The observed 
statistical uncertainty in data is indicated by the arrow, (c) Same as 
(a), but for the /i+jets channel, (d) Same as (b), but for the p+jets 
channel. 

certainties, and the corresponding p value is 0.30. 

We repeat the measurement assuming Gf.s = 1, he-, by eval¬ 
uating the one-dimensional likelihood ±jfs = 1 ), 

and find mf ,ES=1 = 176.88 ±0.41 GeV, which is consistent 
with the final results obtained by maximizing 22 ?(X 7 v;/«f,A:jEs) 
as a function of (/w,,£jes), considering the measured &jes- 
Separating the statistical uncertainty <7,„, into two compo¬ 
nents, one from m, and the other from £jes, we obtain m, = 
174.98 ± 0.41 (stat) ± 0.41 (JES) GeV. 

In contrast to the previous measurement fl4l . we do not use 
the JES determined in exclusive y±jet and dijet events with an 
uncertainty of ss 2% to constrain Ges- We follow this strategy 
because the statistical uncertainty on the measured Gr.s value 
is smaller than the typical uncertainty on the JES and because 
kjEs relates jet energies at detector level to parton energies, 
while JES relates jet energies at detector level to jet energies 
at particle level. 


C. Comparison of data with simulations using the measured 

/, m,, and £jes values 


In this Section, we demonstrate that our simulations agree 
with the data for the measured m, and Gi:s parameters. We use 
a sample of simulated tt events at mf en = 175 GeV and scale 
the energies of jets by Grs = 1-025, while adjusting the mag¬ 
nitude of the three-momentum \p\ to preserve the invariant 
mass of the jet. The fraction of signal events, /, is set to our 
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FIG. 26: The invariant mass of the tt system in the (a) e+jets and (b) /r+jets final states, assuming mf = 172.5 GeV, kj^s = 1, and 
a,j = 7.24 pb~* g0). The same observable in the (c) e+jets and (d) fl +jets final states for mf en = 175 GeV, with energies of jets scaled up by 
fcjES = 1.025, corresponding to their measured values. The fractions of signal events, /, are set to their measured values of 63% in the e+jets 
and 70% in the fl+ jets final states in (c) and (d) (cf. Section [VT^J. The category “Other bgs” encompasses Z+jets, WW, WZ, ZZ, and single 
top quark production, as well as dileptonic tt decays. The production of a IT boson in association with at least one b or c quark is denoted 
as “VT + hf’, while “W + If’ stands for all other flavor configurations. The last bin includes overflow events. The ratio shows the number of 
observed events divided by the number of expected events in a given bin. The band of systematic uncertainty is indicated by the shaded area 
in the ratio plots, and includes contributions from dominant sources: JES, jet energy resolution, jet identification, flavor dependence of the 
detector response, b tagging, lepton identification, lepton momentum scale, trigger, as well as hadronization and higher-order effects for tt 
events. 


measured values of 63% in the e+jets and 70% in the /r+jets 
final states (cf. Section [VI B| i. We consider the same system¬ 
atic uncertainties as in Section [Tv] JES, jet energy resolution, 
jet identification, flavor dependence of the detector response, 
b tagging, lepton identification, lepton momentum scale, trig¬ 
ger, as well as hadronization and higher-order effects for tt 
events. 


Distributions in the invariant mass of the tt system and in 
the invariant mass of the jet pair matched to one of the W 
bosons are shown in Figs. [26] (c-d) and 27 (c-d), respectively, 
together with the predictions from our signal and background 
models. The observables m t j and Mw are reconstructed with 


the same algorithm as described in Section IV As expected, 
the agreement between data and simulations improves com¬ 
pared to Figs. 26 (a-b) and 27 (a-b) when we use the measured 
values of /, m t , and kjES- 


In Fig. [28] we study the distribution in /// , which is de¬ 
fined as the scalar sum of Ej of the four jets, the transverse 
momentum of the tt system, and the transverse momentum 
of the top quarks. The observable p u r is given by the trans¬ 
verse coordinates of the four-momentum vector used to calcu¬ 
late m T f, while is obtained from the four-momenta of final 
state particles used to reconstruct m,-,. These are combined to 
give the four-momenta of the t and t quarks that minimize the 
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FIG. 27: Same as Fig. 26 but for the invariant mass of the dijet system matched to the W boson, assuming mf = 172.5 GeV, Ajes = 1 


and a,j = 7.24 pb 1 |40| in the (a) e+jets and (b) /r+jets final states, and for mf™ = 175 GeV, fejgs = 1.025, and / corresponding to their 
measured values in the (c) e+jets and (d) p + jets final states. 


| m t — nif | mass difference. 

The transverse momenta of the electrons and muons, the 
distributions of rfj, and the transverse mass of the W -+ £V( 


decay are shown in Fig. 29 The invariant of the trijet system 


matched to the hadronic decay of the top quark, reconstructed 
from the same four-momenta used in the reconstruction of m,j 


and My/, is shown in Fig. 30 Notable improvement due to 
reduced backgrounds is observed relative to the Run I mea¬ 
surement I 


Good agreement between data and simulations is observed 
in all variables except ///, where ALPGEN+PYTHIA pre¬ 
dicts fewer events with p n r < 10 GeV than observed in 
data , for whic h a systematic uncertainty is evaluated in Sec¬ 
tion VIIIA 2 The variable /? 7 ° p , for which agreement within 
uncertainties is observed, is shown at detector level. How¬ 
ever, some discrepancy between measured do+/d/?y 0p cor¬ 
rected for detector effects, and SM calculations is present in 
our data I 
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FIG. 28: The scalar sum of Ex of the four jets Hj in (a) e+jets and (b) /l+jets final states, for the transverse momentum of the tt system 
in (c) e+jets and (d) /i+jets final states, as well as for the transverse momentum of the top quarks in (a) e+jets or (b) /i+jets final states. 
All panels are for mf en = 175 GeV, with energies of jets scaled up by Icjes = 1.025, corresponding to their measured values. The fractions of 
signal events, /, are set to their measured values of 63% in the e+jets and 70% in the /i+jets final states in (c) and (d) (cf. Section [VIB| l. The 
list of included background contributions and sources of systematic uncertainties considered are identical to those from Fig. [26] Panels (e) and 
(f) have two entries per event, one for each of the t and t quarks. 
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FIG. 29: Same as Fig.[28] but for the transverse momentum of the (a) electron or (b) muon, for the missing transverse momentum in (c) e+jets 
and (d) fi +jets final states, as well as for the transverse mass of the W —► £V( decay in (e) e+jets and (f) /i+jets final states. 
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FIG. 30: Same as Fig. [28] but for the invariant mass of the trijet 
system matched to the hadronic decay of the top quark in f+jets 
final states. 


VIII. SYSTEMATIC UNCERTAINTIES 


Systematic uncertainties are evaluated using sets of PEs 
constructed from simulated signal and background events, for 
three categories of sources: 

(i) modeling of signal and background events, 

(ii) simulation of the detector response, and 

(iii) analysis procedures and assumptions. 

The contributions from these sources listed in Tableware eval¬ 
uated individually and combined in quadrature to obtain the 
total systematic uncertainty. 

All sources of systematic uncertainties are considered fully 
correlated between e+jets and ft 4-jets final states, except for 
the uncertainties from the lepton momentum scale and the cal¬ 
ibration of the method. In comparison to the previous anal¬ 
ysis using 3.6 fb _1 of integrated luminosity M, the num¬ 
ber of simulated tt events used for the calibration of the ME 
technique and for the evaluation of systematic uncertainties 
has increased significantly since the computation time for Psig 
has been reduced by two orders of magnitude, as discussed 
in Section [VF] As a result, the statistical components of sys¬ 
tematic uncertainties are ss 0.05 GeV in each of the first four 
sources listed in Table[5] namely those from higher-order cor¬ 
rections, initial and/or final-state radiation, hadronization and 
the underlying event (UE), and color reconnection, and are 
« 0.01 GeV for the heavy-flavor scale factor, residual jet en¬ 
ergy scale, flavor-dependent response to jets, lepton momen¬ 
tum scale, jet energy resolution, jet identification efficiency, 
and modeling of multijet events. The statistical components 


are negligible for all the other sources of systematic uncer¬ 
tainty. The uncertainties from the first four sources are evalu¬ 
ated for mf en = 172.5 GeV and = 1 by comparing results 
for m t using different models for signal. All other system¬ 
atic uncertainties are evaluated by applying a recalibration us¬ 
ing simulations that reflect an alternative model to data. The 
statistical components due to the finite number of simulated 
events are never larger than the net difference between the 
central and alternative models for any of the sources of sys¬ 
tematic uncertainty. 

Unless stated otherwise, the systematic uncertainties are 
evaluated using fully simulated events. Background events 
are always included, as described in Sections [VI A| and |VII A| 
The signal fraction used to compose PEs is kept fixed at the 
value obtained for that data-taking epoch and final state, as 
summarized in Table [4j except in evaluating the systematic 
uncertainty from /. The systematic uncertainties from cate¬ 
gory (iii), which are associated with analysis procedures and 
assumptions, are evaluated including the MJ background in 
the ensemble tests. The estimate of all other systematic un¬ 
certainties is not expected to be affected by the inclusion of 
the MJ background contribution in the ensemble tests, which 
is therefore neglected. We quote signed uncertainties for one¬ 
sided sources of uncertainty or sources dominated by a one¬ 
sided component in Table [5] indicating the direction of m t 
change when using an alternative instead of the central model. 
Two-sided uncertainties or uncertainties dominated by a two- 
sided component are marked with a ± or + symbol. For one¬ 
sided sources of systematic uncertainty, the full magnitude of 
the effect is taken in each direction when calculating the total 
quadrature sum to obtain the total systematic uncertainty, in 
accordance with Ref. Q. 


A. Signal and background modeling 


1. Higher-order corrections 


We calibrate the response of the ME technique using tt 
events generated with the LO ALPGEN program, which ac¬ 
counts for real corrections to the Born level process up to 
second-order in a s through the ME, i.e., for 0, 1, or 2 ad¬ 
ditional light partons. To evaluate the effect of higher-order 
virtual corrections on m t , we use signal events generated with 
the MC@NLO program ED, which calculates the processes 
qq —> tt, gg —> tt, and qg —> tt at NLO, as an alternative model. 
Because MC@NLO is interfaced to HERWIG lf63l for simulat¬ 
ing the contributions from leading-logarithmic (LL) correc¬ 
tions and fragmentation to the hard process of interest, we use 
ALPGEN+HERWIG events as the central model for this com¬ 
parison to avoid double-counting the uncertainty from a dif¬ 
ferent showering model. Performing the ensemble tests for 
m t = 172.5 GeV, we find a shift of +0.15 GeV, which we 
assign as the systematic uncertainty due to higher-order cor¬ 
rections. 
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TABLE 5: Summary of uncertainties on the measured m,. One¬ 
sided sources of uncertainty or sources dominated by a one-sided 
component are signed, indicating the direction of m, change when 
using an alternative instead of the central model. Two-sided un¬ 
certainties or uncertainties dominated by a two-sided component are 
marked with a =t or + symbol. 


Source of uncertainty 

Effect on m t (GeV) 

Signal and background modeling: 


Higher-order corrections 

+0.15 

Initial/final state radiation 

+0.06 

Transverse momentum of the tt system 

-0.07 

Hadronization and underlying event 

+0.26 

Color reconnection 

+0.10 

Multiple pp interactions 

-0.06 

Heavy-flavor scale factor 

+0.06 

Modeling of b-quark jet 

+0.09 

Parton distribution functions 

±0.11 

Detector modeling: 


Residual jet energy scale 

±0.21 

Flavor-dependent response to jets 

+0.16 

Tagging of b jets 

+0.10 

Trigger 

±0.01 

Lepton momentum scale 

±0.01 

Jet energy resolution 

±0.07 

Jet identification efficiency 

-0.01 

Method: 


Modeling of multijet events 

+0.04 

Signal fraction 

±0.08 

MC calibration 

±0.07 

Total systematic uncertainty 

±0.49 

Statistical uncertainty 

±0.58 

Total uncertainty 

±0.76 


2. Initial/final state radiation 


The modeling of extra jets from initial/final state radiation 
(ISR/FSR) can affect the measurement of in,. To evaluate the 
contribution from this source, we adjust the amount of radi¬ 
ation by changing the ktfac parameter |64)| that defines the 
renormalization scale in the CKKW scale-setting procedure, 
from its standard value of ktfac = 1 in our ALPGEN+PYTHIA 
simulations described in Section|IIIJ as suggested in Ref. (64) . 
We constrain changes in the amount of ISR/FSR by studying 
the cross section for Drell-Yan events (qq -A Z/y* —» £ + £~) 
measured differentially in the (j)* variable (65), which is re¬ 
lated to the pr of the Z/ y* boson, and has sensitivity to 
the amount of QCD radiation. Our studies, documented in 
Appendix |A) indicate that changes in the ktfac parameter 
by a factor of ±1.5 cover the excursion of MC simulations 
from data, which results in an uncertainty of +0.06 GeV on 
m t . Previously, the systematic uncertainty from modeling of 
ISR/FSR was evaluated using a pure LO generator, PYTHIA, 
and changes made in shower-evolution parameters. The new 
procedure is more accurate, as it uses a larger data sample 
for deriving the ISR/FSR dependence, and employs the same 
generator setup as for simulating the default tt sample. 


3. Transverse momentum of the tt system 


The distribution in the transverse momentum of the tt sys 


tem, pj , is not modeled well as shown in Fig. 28 and a dedi¬ 
cated systematic uncertainty is assigned. This mismodeling 
is corrected by reweighting the contribution from tt events 
differentially in p n r to achieve agreement between the simu¬ 
lations and data. Propagating this effect to m, using ensemble 
tests, we find a shift of —0.07 GeV. 


4. Hadronization and underlying event 


The choice of a model for the parton-shower (PS) evolu¬ 
tion, the hadronization, and the underlying event (UE) can 
lead to a bias on the measured m t value. To evaluate the size 
of this effect, we compare MC samples where PS evolution, 
hadronization, and UE are simulated using PYTHIA with HER- 
WIG as an alternative model. In both cases, the ALPGEN gen¬ 
erator is used to simulate the hard process. The comparison of 
ALPGEN+HERWIG with ALPGEN+PYTHIA will give a convo¬ 
lution of four terms: 

(i) a statistical component due to the limited size of the MC 
samples, 

(ii) a difference in kinematic spectra, most notably in /?", 
between ALPGEN+HERWIG and ALPGEN+PYTHIA, 

(iii) a difference in the response of the detector to jets simu¬ 
lated with HERWIG and PYTHIA, and 

(iv) a difference in the modeling of the PS evolution, 
hadronization, and UE when using ALPGEN+HERWIG 
and ALPGEN+PYTHIA. 

We are interested in the effect from source (iv), since (ii) is al¬ 
ready accounted for in the systematic uncertainty from mod¬ 
eling of ISR/FSR, and (iii) is taken into account in the sys¬ 
tematic uncertainties from modeling of the detector response. 
We reduce the contributions from sources (i), (ii), and (iii) as 
follows: 


(i) The contribution from source (i) is reduced to 
ss 0.05 GeV, about a factor of five smaller than the pre¬ 
vious m, measurement lH4l . by increasing the size of 
MC samples by a factor of « 25. 


(ii) Reshuffling of momenta between the hard and soft 
part of the event (66l causes the p'' r distribu¬ 
tion in ALPGEN+HERWIG to differ significantly from 

To isolate 


ALPGEN+PYTHIA, as shown in Fig. 31 


this effect in the evaluation of the systematic uncer¬ 
tainty, we reweight ALPGEN+PYTHIA in p l j. to match 
the ALPGEN+HERWIG simulation, which results in a 
shift of 0.05 GeV in m,. The uncertainty from mod¬ 
eling of transverse momentum of the tt system already 
accounts for the potential impact on m, of this differ¬ 
ence, as discussed in Section[VIII A 3| 
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FIG. 31: (a) The ratio of predicted differential cross sections from 
ALPGEN+HERWIG to those from ALPGEN+PYTHIA as a function 
of p'j, at particle level. Also shown is the fit with a second-order 
polynomial used for reweighting ALPGEN+PYTHIA events in pj for 
35 GeV < p'j < 200 GeV. Below that range the reweighting is per¬ 
formed according to the histogram, and above with a constant factor 
of 1.67. Uncertainties are only statistical, (b) Same ratio in the range 
p‘j- < 35 GeV. 


and UE by comparing tt events simulated with ALPGEN and 
interfaced to either PYTHIA with DO tune A (the standard MC 
setup), or PYTHIA with the Perugia 2011C tune |67l . While 
the former is based on a Q 2 ordered PS model, the latter is p 2 r 
ordered 12811291 . Furthermore, the latter is a recent tune from 
the Perugia 2011 family that includes LHC data, which is rec¬ 
ommended for this comparison ( 68 ). Both use the CTEQ6L1 
set of PDFs. We verify that the Perugia 2011 family describes 
the transverse jet shapes sufficiently well and is therefore ap¬ 
plicable to Tevatron data dominated by quark jets, despite be¬ 
ing tuned including LHC data dominated by wider gluon jets. 
For simplicity, we carry out the comparison between PYTHIA 
with DO tune A and PYTHIA with Perugia 2011C at detector 
level, and find a change of +0.30 GeV in m t , which serves 
as an upper estimate of the effect evaluated with particle-level 
jet energies, as for the comparison of ALPGEN+HERWIG with 
ALPGEN+PYTHIA. 


(iii) The JES calibration is derived using PYTHIA with 
DO tune A 031 , and is therefore valid only for this spe¬ 
cific configuration. The contribution from source (iii) 
would not be present if a dedicated JES derived with 
HERWIG was available. Like the default JES calibra¬ 
tion using PYTHIA, the alternative calibration using 
HERWIG would be affected by similar uncertainties ac¬ 
counted for in the “detector modeling” category of Ta¬ 
ble [5] Including the difference of the JES derived with 
PYTHIA and HERWIG in the evaluation of the uncer¬ 
tainty from the choice PS evolution, hadronization, and 
UE models would lead to double-counting. To elim¬ 
inate the contribution from the response of the detec¬ 
tor, we use particle-level jet energies. The standard se¬ 
lection as described in Section IV is applied at detec¬ 
tor level to take into account any potential effects from 
the trigger or selection requirements, and the detector- 
level jets are matched to particle-level jets. The match¬ 
ing is done in ( 77 , 0 ) space, using the requirement of 
(./detector? ./particle ) + f^cone /2 = 0.25. The energy of 
each detector-level jet is replaced by the energy of the 
matched particle-level jet, and the magnitude of the 
three-momentum \p\ of the detector-level jet is adjusted 
to preserve its original )], tp, and m parameters. Those 
modified four-momenta of jets are then used, together 
with all other detector-level information in the event, to 
extract m t with the ME technique just as in data. 


Reducing the contribution from source (i), and factoriz¬ 
ing out contributions from sources (ii) and (iii), as described 
above, we find an effect of +0.26 GeV from the choice of the 
model for PS evolution, hadronization, and UE. We verify that 
^jes GEN+HERW 1G anc * k^| GEN+PYTHIA are consistent within uncer¬ 
tainties, as expected at particle level, indicating that the effect 
of JES has been factorized out. 

As a cross-check, we also evaluate the impact of a differ¬ 
ent PS model and a change in the models for hadronization 


5. Color reconnection 


Our standard PYTHIA tune, DO tune A, does not include 
an explicit model for color reconnection (CR). Color recon¬ 
nection is sometimes confused with strong color correlation 
between the hard scattering process and the underlying event. 
The latter is an integral part of PYTHIA DO tune A, as not in¬ 
cluding it would violate confinement. Recently, research has 
been carried out to include CR in PYTHIA l69l . To evalu¬ 
ate the CR effect, we compare identical ALPGEN events inter¬ 
faced to PYTHIA with two different tunes, Perugia 2011 and 
Perugia 2011NOCR. While the former tune includes an ex¬ 
plicit CR model, the latter does not. The Perugia 2011 tune 
rather than Perugia 2011C is used as an alternative model 
because Perugia 2011 uses the same set of CTEQ5L PDFs 
as the reference model Perugia 2011NOCR, while Perugia 
2011C uses CTEQ6L1. This choice of generator setup is 
advised in Ref. ( 68 l . since Perugia includes a better, more 
phenomenologically-motivated CR model than the ACRpro 
tune ED- The new tune models the color-string survival prob¬ 
ability depending on the distance in rapidity between the be¬ 
ginning and the end of the color string, whereas the previous 
model uses only an overall approximation. We find an effect 
of +0.10 GeV on m t from this source of systematic uncer¬ 
tainty. 


To cross-check the impact of modeling of multiple parton 
interactions (MPI) on the measurement, we compare events 
generated with ALPGEN interfaced to PYTHIA with the Pe¬ 
rugia 2011mpiHi and Perugia 2011 tunes, where the former 
features harder MPI, albeit at a smaller rate, compared to the 
latter. We find an effect of 0.14 GeV, with a statistical compo¬ 
nent of « 0.15 GeV due to limited MC event statistics. Since 
HERWIG and PYTHIA use different models for MPI, an uncer¬ 
tainty from the choice of a model for MPI is already included 
in the uncertainty from modeling of hadronization and the UE 
discussed in Section VIIIA 4| In addition, it is at least partly 
included in the comparison between the Perugia 2011 and 
Perugia 2011NOCR tunes discussed above, since the latter. 
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without CR, must simultaneously adjust the amount of MPI 
to describe data. Therefore, the comparison between Perugia 
201 ImpiHi and Perugia 2011 is treated as a cross-check. 


6. Multiple pp interactions 

Multiple pp interactions can potentially influence the mea¬ 
surement of m t . We reweight the spectrum in instantaneous 
luminosity of our simulated MC samples to the instantaneous 
luminosity profile found in data for the respective data-taking 
epoch and final state. This default correction is derived for 
each data-taking stream, while the instantaneously luminos¬ 
ity profile in a subset of this data can differ due to the trigger 
and selection requirements. To evaluate this systematic uncer¬ 
tainty, we obtain an alternative calibration where an additional 
reweighting in instantaneous luminosity is applied after all se¬ 
lections to bring its spectra in tt simulations and in data into 
agreement. We find a change in m, of —0.06 GeV. 


7. Heavy-flavor scale factor 

As mentioned in Section [IIIJ a heavy-flavor scale factor of 
kbb,cc = 1-47 is applied to the Wbb and Wcc cross sections to 
correct the heavy-flavor content of the W+jets MC samples to 
expectations from NLO calculations. The uncertainty on kbh, a 
is « 15% |34) . Moreover, there is an additional uncertainty of 
ss 12% on the W + c content of the sample [34). These uncer¬ 
tainties are combined linearly and rounded to 30%. The effect 
°f kbh.cc variations is propagated to m, using PEs. This in¬ 
troduces variations in the kinematic distributions of the com¬ 
bined contribution from W +jets production. As we are in¬ 
terested only in the effect of the changed background com¬ 
position, we do not change / when constructing the pseudo¬ 
experiments. We find an effect of +0.06 GeV due to this 
source of systematic uncertainty. 


8. Modeling ofb-quark jet 

Uncertainties in the simulation of b quark fragmentation 
can affect the measurement of m t through several aspects 
of the analysis such as b tagging and the TF. Such effects 
are studied by reweighting the simulated tt events according 
to other possible fragmentation models for b quarks. The 
standard tt MC samples consist of events reweighted from 
the standard PYTHIA b fragmentation function to a Bowler 
scheme 1711 that has been tuned to LEP (ALEPH, OPAL, and 
DELPHI) data |[72l . To evaluate the potential bias from this 
source, tt events are further reweighted to account for dif¬ 
ferences in SLD |72) and LEP data, resulting in a shift of 
+0.08 GeV in the measured m t . 

The detector response to ft-quark jets can be different in 
the presence of semileptonic decays of b or c quarks. The 


incorrect simulation of semileptonic branching ratios in b and 
c quark decays can thus result in a systematic effect on m t . We 
take an uncertainty of 0.05 GeV determined in Ref. Il48l as the 
contribution from this source. Quadratically combining the 
two sources of uncertainty, we find a shift in m t of +0.09 GeV. 


9. Part on distribution functions 

To evaluate the uncertainty from the choice of PDFs, the 
standard ALPGEN+PYTHIA sample at mf en = 172.5 GeV is 
reweighted to match possible excursions in the PDF parame¬ 
ters provided in CTEQ6M H33I . Ensemble tests are repeated 
for each of these changes and the total uncertainty is evaluated 
using the following formula li33) : 

1 f 20 

n£ 

where the sum runs over PDF uncertainties in the positive 
(S+) and negative (S- ) excursions. We find changes in m, 
of ±0.11 GeV due to this source of systematic uncertainty. 
The effect on m, from replacing the central PDF set from 
CTEQ6M by that from CTEQ6L1 is much smaller. 


B. Detector modeling 

1. Residual jet energy scale 

The overall jet energy scale factor &jes is fitted simultane¬ 
ously with m t , and the corresponding contribution to the sta¬ 
tistical uncertainty is included in a ,„ t . However, Lips does not 
account for possible effects from uncertainties in the JES cor¬ 
rections that are differential in E or /] of the jet. To estimate 
the potential bias on m t from the residual JES uncertainty, we 
use a variety of approaches: 

(i) As for the 3.6 fb 1 analysis, we parameterize the de¬ 
pendence of the upper side of the uncertainty of the JES 
on energy E of the jet, in four bins in rj, as in Ref. 02). 
This parameterization is applied to scale up jet energies 
and to provide an alternative calibration. A single con¬ 
stant offset is added to the parameterization across all 
four rj bins, such that the average correction applied to 
the energies of all the jets in the sample is zero to avoid 
double-counting of acceptance effects. We find an ef¬ 
fect of +0.13 GeV on m t . 

(ii) We also apply the downward changes, and find an effect 
of —0.08 GeV on m t . 

(iii) We apply the positive JES changes directly, i.e., the 
momentum of each jet is increased by its 1 SD uncer¬ 
tainty, without any parametrization. In contrast to (i) 
and (ii), the jet energies are not rescaled on average. 


-m2 


[Am t (SJ)-Am t (S i )] 


(38) 
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and the method is expected to account for the change in 
the overall JES through the £jes parameter in the two- 
dimensional likelihood fit, at the cost of some double¬ 
counting of acceptance effects. We find a change in m, 
of —0.20 GeV with this procedure, while £jes changes 
by —0.018, as expected. 

(iv) The energies of jets are increased using a linear param¬ 
eterization, which is unity for E = 0 GeV and increases 
linearly with E, such that it tangentially approaches the 
upper limit of uncertainty of the JES. The linear param¬ 
eterization is obtained separately for each data-taking 
epoch, and changes slightly due to changes in the JES 
uncertainty. No overall rescaling of the jet energies is 
applied. We find —0.21 GeV with this approach for the 
impact on m t . 

As only one dependence of the true JES on E and rj of the jet 
is realized within the uncertainty corridor, we take the max¬ 
imum envelope of the above uncertainties of 0.21 GeV and 
assign it as the systematic uncertainty due to the residual JES 
correction. In previous analyses MM, only approach (i) 
was used, which corresponds to an uncertainty of 0.13 GeV 
using our results. 


2. Flavor-dependent response to jets 


The standard JES calibration used in this measurement con¬ 
tains a flavor-dependent jet response correction, which is re¬ 
sponsible for bringing the simulation of calorimeter response 
to jets into agreement with that observed in data for a given 
jet flavor m. The flavor-dependent correction is given by 


1 L iEi-Rf** 

~ (F)r+jet' L.Ei-R^’ 


(39) 


where the index i runs over all constituent particles that be¬ 
long to a particle jet, and Rj is the detector response to particle 
i, given differentially in pj and /] of that particle. The factor 
(F) r+jet ensures that the flavor-averaged part of the JES cali¬ 
bration, extracted from exclusive 7 +jet events, is preserved. 
The correction F con is largest for +quark jets and gluon jets, 
and is relatively small for light quark jets, compared to the 
flavor-averaged JES, as shown in Fig. 34 of Ref. 03. We are 
in a position to apply such a correction in a straightforward 
way, since the JES is corrected to particle level both for data 
and MC. To propagate the effect of the uncertainty on F corr to 
our measured m, value, we change F con - by ±1 SD for each 
data-taking epoch and final state. We apply this factor to the 
jets in our signal MC samples to get alternative calibrations, 
resulting in a systematic uncertainty of +0.16 GeV after sym- 
metrization. This uncertainty accounts for the difference in 
detector response to different jet flavors, in particular fc-quark 
jets versus light quark jets. 

We perform additional studies beyond those done in 
Ref. fl5ll to verify that F corr is correctly estimated: 


(i) The derivation of F corr requires a priori input for the 
relative contributions of b quark, gluon, and light quark 
jets to the samples of y+jets and dijet events that are 
used for its derivation. A variation on the assumption 
of these relative contributions would result in somewhat 
different fitted RJ lala in Eq. (391. The relative contribu¬ 


tions from different jet flavors to the sample are taken 
from PYTHIA with DO tune A. To estimate the effect of 
this assumption, we exploit the fact that y+jets and di¬ 
jet samples show a rather different sample composition 
in terms of jet flavor, and obtain F corr in those two sam¬ 
ples independently, rather than simultaneously in the 
standard procedure. We find that the results for F C0IT in 
y+jets and dijet events are consistent with each other 
within assigned uncertainties. 

(ii) Similarly, a priori input about the particle composition 
of the jet is needed for the derivation of F COIT , which is 
also taken from PYTHIA with DO tune A. A variation of 
the composition of the jet would result in a change to 
the numerator and the denominator of Eq. because 
the sum in i runs over all the particles in the jet, and 
to Rf ata , which are determined from a fit to data in y + 
jets and dijet events. To estimate an upper limit on the 
impact of the particle composition on F con - for />quark 
jets, for which F con - is largest, we use simulated tt events 
to study the ratio 


0 = 
^corr 


x COIT, ALPGEN+HERWIG 

F b 

^COIT, ALPGEN+PYTHIA 


(40) 


where F* )rr alpgen+P ythia is the default correction from 

Ref. (15), and F C orr, alpgen+herwig takes the particle com¬ 
position of the jet predicted by ALPGEN+HERWIG as in¬ 
put, but uses the same Rf ata and R^ ]( as in the default 
case. The ratio Q b on is shown for f+jets final states in 


Fig. 32 The excursion of the ratio from unity is within 


0.5% for the pr spectrum of /i-quark jets in tt events. 

The ratio Q b orr can only provide an upper limit on the 
impact of the assumption of the particle composition 
of the jet on the flavor-dependent calorimeter response, 
as it does not involve a coherent JES calibration using 
y+jet and dijet events. In particular, R^ 213 , which are 
determined with PYTHIA, are unchanged, while a cali¬ 
bration of the single particle responses using HERWIG 
would bring F C orr, alpgen+herwig into better agreement 
with y+jet and dijet data, which F b orr alpgen+pythia is 
also tuned to, resulting in a Q b ort that is closer to unity. 
Given these findings, we conclude that the uncertainty 
on m t from flavor dependence of the calorimeter re¬ 
sponse to jets of +0.16 GeV covers the effect from 
different parton shower and hadronisation models in 
PYTHIA and HERWIG. Taking our upper limit estimate 


using Eq. (40 1 at face value would increase the system¬ 
atic uncertainty on m, from flavor dependence of the 
calorimeter response to jets to +0.24 GeV, and the to¬ 
tal uncertainty of the measurement to 0.78 GeV. 
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(iii) We perform another cross check of the energy scale for 
/;-quark jets using a similar approach to Ref. fim by 
studying the ratio of transverse momenta of /^-tagged 
jets to jets without a b tag. For events with > 2 ^-tagged 
jets, this ratio is defined as 


Rbq = 


Pt +Pt 
Pt +P? ’ 


(41) 


where bj refers to the ^-tagged jets with the highest pr, 
and similarly q, for the other two jets. In events with 
one i>-tagged jet, the ratio 


Rbq = 


2 p b T 


Pt +Pt 2 


(42) 


is computed. The distribution of Rq q after all selections 


is shown in Fig. 33 for m t = 175 GeV and a t j = 7.7 pb. 


as provided by the ME technique. Good agreement be¬ 
tween data and simulations is observed. We extract the 
overall b quark JES kbJES by performing a template fit 
to the distribution in R/ Mr and obtain 


JfcbJES = 1 -008 ± 0.0195 (stat) (syst). (43) 

The systematic uncertainties include dominant sources 
like JES, jet energy resolution, jet identification, fla¬ 
vor dependence of the detector response, and b tagging, 
which provides the main contribution. The measured 
kbjES value is consistent with unity, and underlines the 
correctness of F colr from Eq. 39 


3. Tagging of b jets 


Discrepancies in the b tagging efficiency between data and 
MC simulations can lead to a systematic shift in the extracted 
m t . To evaluate the effect of possible discrepancies, the stan¬ 
dard corrections to b tagging rates of ii.d.s quark jets and glu¬ 
ons, c quark jets, and fc-quark jets, which depend on rj and 
Pt, are changed within their uncertainties. We evaluate the 
uncertainty from those corrections by simultaneously decreas¬ 
ing the b tagging efficiencies for b and c quark jets within their 
uncertainties, while increasing the efficiencies for jets from 
all other quarks, and find +0.10 GeV after symmetrizing their 
impact on m,. 


t 1.01 

o 

o 

1.005 

1 

0.995 

0.99 


D0 MC + 

l+jets 


Simulation 
Uncertainty 
Ideal RL„ 



50 100 150 

pf [GeV] 


FIG. 32: The ratio Ocorr — ^corr, alpgen+herwigAcoit.alpgen+pythia 
for £>-quark jets in simulated tt events and f+jets final state as a 
function of the pj of the jet. The ideal <2corr °f unity is indicated 
as the broken line. The uncertainties on the points reflect the finite 
size of the MC sample. The total uncertainty on the standard flavor- 
dependent correction for ALPGEN+PYTHIA is indicated as the band 
around unity. 


5. Lepton momentum scale 


The momentum scale of electrons and muons is known with 
some finite precision. To evaluate the impact of this effect, we 
change the momentum scale of a given lepton flavor, while 
keeping the momentum scale of the other flavor constant, and 
combine in quadrature the effects found for the two flavors. To 
evaluate the impact for a given lepton flavor, we assume the 
worst case, where the momentum scale is shifted in the same 
direction for all data-taking epochs. We perform this for the 
up and down changes, and find an effect that is consistent with 
zero within uncertainties for electrons, and corresponding to 
0.01 GeV for muons. Combining the two values in quadrature 
yields 0.01 GeV, which we quote for this source of systematic 
uncertainty. 


6. Jet energy resolution 


4. Trigger 

To evaluate the impact of the trigger on our analysis, we 
apply trigger weights that simulate the impact of the lepton 
plus jet trigger differentially in pt, R, and (f) of the lepton. An 
alternative calibration with those trigger weights applied is ob¬ 
tained, resulting in a change in the extracted m t of ±0.01 GeV. 


An additional smearing of jet energies is applied to all MC 
samples as part of the JES calibration in order to achieve bet¬ 
ter agreement of MC simulations and data m. To evalu¬ 
ate the possible effect of disagreement between data and MC 
simulations in jet energy resolutions on m t , we produce sim¬ 
ulated MC samples where the jet energy resolution correction 
is changed up and down by its uncertainty. Event probabili¬ 
ties are then re-calculated and ensemble tests are repeated. We 









31 



FIG. 33: The ratio Rf, q = (p* 1 + P* 2 )/ (Pj +Pj) f° r events with 
> 2 fc-tagged jets, and Rj )q = 2p*/(p^i + pj) for events with one 
^-tagged jet, after all selections. The category “Other bgs” encom¬ 
passes Z+jets, WW, WZ, ZZ, and single top quark production, as 
well as dileptonic tt decays. The production of a IT boson in asso¬ 
ciation with at least one b or c quark is denoted as “Vk + hf”, while 
“W + If’ stands for all other flavor configurations. The last bin in¬ 
cludes overflow events. The ratio shows the number of observed 
events divided by the number of expected events in a given bin. The 
band of systematic uncertainty is indicated as a shaded area in the 
ratio plots, and includes dominant sources: JES, jet energy resolu¬ 
tion, jet identification, flavor dependence of the detector response, 
and b tagging. 


C. Method 


1. Modeling of multijet events 


Since we require four jets, and at least one of them b- 
tagged, the contribution of the MJ background to our sample 
of tt candidate events is small (cf. Table [2}. Nonetheless, we 
explicitly account for the next-to-dominant contribution from 
this background when constructing the PEs for the calibration 


of the ME technique, as described in Sections VIA and VIIA 


This was not done in the previous measurement using 3.6 fb 
of integrated luminosity To evaluate a systematic uncer¬ 
tainty from the modeling of MJ background, we do not in¬ 
clude it in the ensemble tests for an alternative calibration in 
m t , which yields a change in m, of +0.04 GeV. 


2. Signal fraction 


We apply the / values measured in each data-taking epoch 
and final state, as summarized in Table |4j to construct the PEs 
used for calibrating the response of the ME technique in m, 
and kjEs ■ To evaluate the systematic uncertainty from /, we 
use alternative calibrations in m, and £jes, where / is changed 
by ±5%, simultaneously for all data-taking epochs and final 
states. The 5% value is motivated by the magnitude of the 
systematic uncertainty on the measured tt production cross 
section in the DO data ED, ignoring the uncertainty from in¬ 
tegrated luminosity. We obtain ±0.08 GeV as an effect from 
this source of systematic uncertainty. 


3. MC calibration 


find an effect of ±0.07 GeV on m, from jet energy resolution, 
after symmetrizing the excursions. 


7. Jet identification efficiency 

The efficiency for the identification of jets is slightly higher 
in MC simulation than in data, and is therefore corrected as 
part of the JES calibration ED to achieve better agreement 
between MC simulations and data in jet identification effi¬ 
ciencies. To evaluate the potential impact from this source on 
m t , the jet identification efficiencies in the signal MC sample 
are decreased according to their uncertainties. Subsequently, 
event probabilities are recalculated and ensemble tests are re¬ 
peated. The extracted value of m t using this alternative cali¬ 
bration differs from the central value by —0.01 GeV, which is 
quoted as systematic uncertainty. 


The statistical uncertainties due to the limited size of MC 
samples that are used to construct PEs to study the response 
of the ME technique, as discussed in Sections fVI A| and |VII A| 
contribute an uncertainty on the final calibrated value of m, 
that is statistical in nature. To determine this contribution, 
we note the uncertainties on the offset and slope parameters 
of the linear calibrations in Figs. im and [12] and propagate 
them, assuming Gaussian uncorrelated uncertainties, to the 
extracted m, values in each of the data-taking epochs and fi¬ 
nal states. Combining the individual m t results for the up and 
down changes and symmetrizing, we find an uncertainty of 
±0.07 GeV due to the limited size of MC samples. 


D. Summary of uncertainties 

Quadratically combining contributions from all sources of 
systematic uncertainties from Table [5] we obtain a total sys¬ 
tematic uncertainty of 0.49 GeV. The total uncertainty of 
the measurement of 0.76 GeV is obtained by combining the 
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July 2014 (‘preliminary) 

CDF-II dilepton * 

170.80 ±3.26 (±1.83 ±2.69) 

D0-II dilepton 

174.00±2.80 (±2.36+ 1.49) 

CDF-II alljets * 

175.07±1.95 (±1.55+ 1.19) 

CDF-II MET+Jets 

173.93±1.85 (±1 26+ 1.36) 

CDF-II lepton+jets 

172.85±1.12 (±0.52+0.98) 

D0-II lepton+jets 

174.98 ±0.76 (±0.41 + 0.63) 

Tevatron combination * 174 . 34 ±o,64 (± 0 . 37 + 0 . 52 ) 

(Run 1 and Run II) (±stat ± syst) 


% 2 /dof = 10.8/11 (46%) 

_1_ 1 _ 


165 170 175 180 185 

m t [GeV] 


FIG. 34: Summary of the measurements performed in Run II of 
the Tevatron which are used as inputs to the Tevatron combina¬ 
tion GO- The inner error bars indicate the statistical uncertainties, 
while the outer bars represent the total uncertainties. The Tevatron 
average value of m t obtained using input measurements from Run I 
and Run II is given at the bottom and its uncertainty is shown by the 
band. The Figure is adapted from Ref. (8). 


total systematic uncertainty and the statistical uncertainty of 
0.58 GeV. The statistical uncertainty consists of two parts: 
0.41 GeV fromm, alone and 0.41 GeV from Gi:s- Considering 
the latter as a source of systematic uncertainty, as is done for 
the Tevatron combination |8j. the total systematic uncertainty 
including the statistical contribution from the in situ constraint 
on the JES is 0.64 GeV. 


IX. COMPARISON WITH PREVIOUS MEASUREMENTS 

A. Comparison with the Tevatron average 

Our result is included in the Tevatron combination from 
July 2014 |8|, which takes into account 10 published and 2 
preliminary results from the CDF and DO collaborations us¬ 
ing pp collision data from Run I and Run II of the Fermilab 
Tevatron Collider. Considering correlations between sources 
of systematic uncertainty as described in detail in Ref. (7J, the 


final result is m t = 174.34 ± 0.64 GeV, with a relative weight 
of 67% from the DO measurement. An overview of the in¬ 
put measurements performed using Run II data is presented in 
Fig. [34] All measurements from Run I and Run II are consis¬ 
tent with x 2 = 10.8 for 11 degrees of freedom, which corre¬ 
sponds to a x 2 probability of 46%. 


B. Comparison with the world average 

Our result can be compared with the current world aver¬ 
age of mt = 173.34 ±0.76 GeV 0 , which encompasses 11 
measurements from the ATFAS, CDF, CMS, and DO Collab¬ 
orations, excluding this measurement. We do this compari¬ 
son considering the full uncertainty in the world combination 
and only the statistical uncertainty in our measurement, which 
provides a reasonable approximation, given the large correla¬ 
tion among the experiments for most sources of systematic 
uncertainty. Furthermore, we use the simplifying assumption 
that there is no correlation between the statistical uncertainty 
of our new measurement and of the world average. We find 
consistency at the 1.7 SD level. Due to the complicated corre¬ 
lation of systematic uncertainties, a more detailed comparison 
should be performed in a separate document, with participa¬ 
tion from all the collaborations supplying measured m t val¬ 
ues as inputs, and include all updates since the publication of 
Ref. 0. 


X. SUMMARY 

In summary, we have performed a measurement of the mass 
of the top quark using the matrix element technique in tt 
candidate events in lepton±jets final states using 9.7 fb 1 of 
Run II integrated luminosity collected by the DO detector at 
the Fermilab Tevatron pp Collider. The result, 

m t = 174.98 ±0.58 (stat + JES) ±0.49 (syst) GeV, or 
m t = 174.98±0.76 GeV, 

is consistent with the values given by the current Tevatron and 
world combinations of the top-quark mass mm and achieves 
by itself a similar precision. With an uncertainty of 0.43%, 
it constitutes the most precise single measurement of the top- 
quark mass, and is « 70% more precise than the next-to-most 
precise single measurement 174). The total systematic uncer¬ 
tainty of our result is smaller than that of any other single 
measurement. 
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Appendix A: Constraining the amount of radiation using 
Drell-Yan events 


Because of the limited statistics in tt candidate events 
recorded by DO in Run II, it is difficult to use tt events di¬ 
rectly to study ISR and FSR. We therefore take an alternative 
approach using the Z/y* —> £ + £~ process to experimentally 
constrain ISR and FSR. In the following, we motivate this ap¬ 
proach, briefly review the measurement of Z/y* —> £ + £ pro¬ 
duction li65l that we use to obtain the dependence of ISR and 
FSR, and present our results. 


1. Motivation 


At the Tevatron, It pairs are dominantly produced via qq 
annihilation. The ISR will stem therefore dominantly from the 
process where an initial-state quark radiates a gluon, which 
can be described by the Pq^> qg splitting function. In the LO 
picture of the tt process, the final-state partons are also quarks, 
and FSR is therefore described by the same Pq^ qg splitting 
function. At the same time, the production of Z bosons at 
the Tevatron is also dominated by qq annihilation, and ISR 
in Z/y* —> £ + £~ can be also described by the Pq^ qg splitting 
function. A measurement of ISR in Z/y* — > t + t~ events can 
therefore be used to constrain ISR and FSR in tt events. 

The experimental approach outlined above was first pro¬ 
posed in Ref. ED. Beyond Ref. 11751 . we use the 0* vari¬ 
able 1651 . which is more sensitive to soft ISR. Since this mea¬ 
surement is corrected for detector effects, it can be used to 
compare directly generated Z/y" —> £ + £~ events with data. 


2. The measurement of the Z/y* —» (.+£- cross section 


The 0* observable used to establish sensitivity to ISR in 
Ref. lf65l is defined as 

0* = tan (0 acop /2) sin(0*), (Al) 

where 0 acO p is the acoplanarity angle, given by 0 acop = Tt — 
Atj) 1 1 , and Ad' 1 is the difference in azimuthal angle be¬ 
tween the two lepton candidates. The variable 0/ is a measure 
of the scattering angle of the leptons relative to the proton 
beam direction in the rest frame of the dilepton system, and is 
defined as 


cos(0*) = tanh ^ ^ ^ , (A2) 

where 7 ] and t ] + are the laboratory pseudorapidities of the 
negatively and positively charged lepton, respectively. 

The measurement of d<T z /y,/dtp* is corrected for de¬ 
tector effects, and performed in the fiducial region defined by 
the presence of the two oppositely charged leptons, as sum¬ 
marized in Table [ 6 ] The electron four-momentum is defined 


TABLE 6 : Summary of kinematic requirements used to define the 
measurement of dcr z /y,_>£+e-/d(j>* in the fiducial region of Ref. 1651. 


Charged leptons: 
£ + t~ system: 


p T > 20 GeV (771 < 1.1 or 1.5 < |tj| < 3.0 (e) 
Pt > 15 GeV j 77 j <2.0 (ju) 

60 GeV < m e+( - < 120 GeV 


as the sum of the four-momenta of all electrons and photons 
in the A R < 0.2 cone around the detector-level electron. No 
cone summation is applied to muons. The measurement is per¬ 
formed in three bins of rapidity of the boson of 0 < bl < 1 . 
1 < |v| < 2, and |y| > 2 in e + e , and in two bins of rapidity, 
0 < |y| < 1 and 1 < |y| < 2 in final states. This en¬ 

sures that Pq->q g is probed in different kinematic regimes of 
the square of the four-momentum transfer Q 2 . 


3. Results 


As described in Sec. [HI] ALPGEN+PYTHIA is used to gen¬ 
erate tt MC events, where Feynman diagrams accounting for 
hard radiation from the initial or final state partons up to 
second-order in a s are introduced through the hard ME in 
ALPGEN, i.e., in bins of 0, 1, and 2 additional light partons, 
which are then combined according to their cross sections. 
For consistency, we use an identical setup for generating Z —I 
£ + £~ events. The ratio of unfolded da z /y*^.(+(- /dtj)* data 
to the predictions from ALPGEN+PYTHIA at particle level, is 
shown in Fig. 35 for the <? + e“ and in Fig. 36 for the 
final states. The ALPGEN+PYTHIA model is not able to de¬ 
scribe the data over the entire spectrum in (j)*, but the overall 
agreement is reasonable. In the analysis, we identify up and 
down changes in the radiation rate that cover the excursions 
of nominal MC predictions from the data. 

Modifying parameters such as Aqcd or a s adjusting the rate 
of ISR/FSR in the PS, as is done in stand-alone PYTHIA, leads 
to unphysical results for ALPGEN+PYTHIA, since the MLM 
matching of the PS to the hard ME calculation partly compen¬ 
sates for such effects in the PS (64). Thus, the preferred ap¬ 
proach is to change the amount of ISR/FSR through rescaling 
of the renormalization scale in the CKKW scale-setting proce¬ 
dure ktf ac from its default of unity, as suggested in Ref. (64) . 
Our studies indicate that changing the ktf ac parameter by a 
factor of ± 1.5 provide coverage for the excursions of the nom¬ 
inal MC predictions from the data. This is demonstrated in 
Fig. 35 for the e + e~ and in Fig. 36 for the final states. 

Based on these studies, we evaluate the systematic uncertainty 
from ISR/FSR in tt events by setting the ktf ac parameter to 
2/3 for the up and 3/2 for the down point. 
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FIG. 36: (a) Ratio of the unfolded differential cross 

section d<T z /y«—/d0 * and the nominal predictions from 
ALPGEN+PYTHIA for 0 < |y| < 1, where the error bars represent 
the total uncertainties in the data. The ratios of predictions obtained 
with ALPGEN+PYTHIA where the amount of radiation is changed up 
and down relative the nominal predictions from ALPGEN+PYTHIA 
are shown, respectively, as red broken and blue dash-dotted lines. 
The statistical uncertainties from the finite number of simulated 
MC events are negligible compared to the uncertainties in the data, 
(b) Same as (a), but for 1 < \y\ < 2. 


FIG. 35: (a) Ratio of the unfolded differential cross 

section d(7 z /y,^ e+e - /d(j>* and the nominal predictions from 
ALPGEN+PYTHIA for 0 < |y| < 1, where the error bars represent the 
total uncertainties in the data. The ratios of predictions determined 
with ALPGEN+PYTHIA where the amount of radiation is changed up 
and down relative the nominal predictions from ALPGEN+PYTHIA 
are shown, respectively, as red broken and blue dash-dotted lines. 
The statistical uncertainties from the finite number of simulated 
MC events are negligible compared to the uncertainties in the data, 
(b) Same as (a), but for 1 < y| < 2. (c) Same as (a), but for |y| > 2. 
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Appendix B: Transfer function parameters 


In this Appendix we summarize the jet transfer function parameters a,- and b, defined in Eq. GD as a function of data-taking 
epoch, jet flavor, and pseudorapidity bin. Jet flavor is denoted as “light jets” for m , of, c,s quark jets, “b^ jets” for /;-quarl< jets with 
a soft muon tag, and “b jets” for all other /;-quark jets. The parameter <33 is fixed to 0 in order to improve the convergence of the 
double-Gaussian fit function. The parameter values are given in units of GeV (a \), 1 (b\), GeV ( 122 ), 1 (bo), 1 (A 3 ), GeV 1 (/?;), 
GeV (o 4 ), 1 (b 4 ), GeV (a 5 ), and 1 (b 5 ). 


TABLE 7: Jet TF parameters a, and b, for Run Ha and |rj| < 0.4. 


i 

Light jets 

cij bj 

b jets 

a, bj 

b^ jets 
di 

b, 

1 

-1.29 x 10 +u 

2.21 x 10~ 2 

2.92 x 10 +u 

-1.91 x 10 _1 

6.06 x 10 +u 

-1.40 x 10 1 

2 

4.04 x 10+° 

1.09 x 10-' 

2.77 x 10+° 

2.07 x lO^ 1 

1.67 x 10 +0 

1.65 x 10~ 1 

3 

0 

2.35 x 10~ 4 

0 

4.09 x 10~ 2 

0 

8.68 x 10- 5 

4 

2.37 x 10 +1 

-1.70 x 10-' 

-9.96 x 10+° 

4.68 x 10~ 2 

4.60 x 10+ 1 

-4.15 x lO^ 1 

5 

1.89 x 10 +1 

9.64 x 10~ 2 

3.82 x 10+° 

9.23 x 10~ 2 

1.83 x 10+ 1 

1.44 x 10 _1 




TABLE 8: Jet TF parameters a; and b, for Run Ha and 0.4 < 1 77 1 

<0.8. 


i 

Light jets 

di bj 

b jets 

a i bi 

bf! jets 

a,- bj 

1 

-6.87 x 10 1 

8.60 x 10~ 4 

4.34 x 10+° 

-2.11 x 10-' 

6.88 x 10 +o 

— 1.56 x 10-' 

2 

3.64 x 10+° 

1.19x10-' 

3.13 x 10+° 

2.03 x 10- 1 

1.71 x 10+° 

1.65 x 10-' 

3 

0 

2.73 x 10~ 4 

0 

3.40 x 10- 2 

0 

9.12 x 10- 5 

4 

2.48 x 10 +1 

-1.73 x 10-' 

-8.77 x 10+° 

2.43 x 10- 2 

4.71 x 10 +1 

-3.41 x 10-' 

5 

1.89 x 10+ 1 

9.84 x IQ- 2 

3.45 x 10+° 

9.94 x IQ- 2 

2.09 x 10 +1 

1.07 x 10-' 


TABLE 9: Jet TF parameters a, and bj for Run Ha and 0.8 < | rj| < 1.6. 

i 

Light jets 

bj 

b jets 

cij b. 

bji jets 

G[ 

bj 

1 

4.53 x 10+° ” 

-7.58 x 10 -2 

1.29 x 10+' 

-1.58 x 10-' 

6.00 x 10 +u 

-1.08 x 10 1 

2 

3.32 x 10+° 

1.55 x 10-' 

5.24 x 10+° 

2.40 x 10- 1 

6.04 x 10+° 

1.30 x 10 1 

3 

0 

4.28 x 10~ 3 

0 

1.76 x 10- 1 

0 

1.00 x 10- 3 

4 

9.02 x 10+° 

7.90 x 10- 2 

-5.98 x 10+° 

-6.28 x 10- 2 

6.00 x 10+ 1 

-7.50 x 10-' 

5 

1.30 x 10+' 

8.39 x 10“ 2 

1.63 x 10+° 

1.82 x IQ- 1 

1.00 x 10 +1 

2.00 x 10- 1 


TABLE 10: Jet TF parameters a, and bj for Run Ha and 1.6 < | rj | < 2.5. 


i 

Light jets 

cij bj 

b jets 

a, bj 

bfj jets 
ai 

bj 

1 

1.62 x 10+ 1 

-2.28 x 10-' 

1.07 x 10+' 

-3.07 x 10- 1 

3.28 x 10+' 

-4.28 x 10-' 

2 

2.78 x 10+° 

1.46 x 10-' 

3.61 x 10+° 

1.23 x 10- 1 

8.34 x 10+° 

7.76 x 10- 2 

3 

0 

6.15 x 10- 3 

0 

6.98 x 10- 3 

0 

7.06 x 10- 3 

4 

1.86 x 10+ 1 

-1.30 x 10- 3 

5.35 x 10+° 

-4.33 x 10- 2 

2.07 x 10 +1 

-1.31 x 10 1 

5 

1.47 x 10+' 

8.60 x IQ- 2 

1.38 x 10+' 

9.66 x 10- 2 

1.18 x 10 +1 

1.08 x 10- 1 
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TABLE 11: Jet TF parameters a,- and b, for Run Ilbl and J r\ \ < 0.4. 


i 

Light jets 

a,- bj 

b jets 

fl/ 

bi 

bn jets 

a, b. 

i 

-2.17 x 10+° 

2.39 x 10~ 2 

8.74 x 10-' 

-1.86 x 10-' 

5.26 x 10+° 

-1.29 x 10-' 

2 

4.08 x 10+° 

6.66 x 10- 2 

2.67 x 10+° 

1.82 x 10- 1 

1.45 x 10+° 

1.51 x 10-' 

3 

0 

2.77 x 10~ 4 

0 

4.41 x 10- 2 

0 

2.13 x 10- 4 

4 

1.27 x 10 +1 

-1.76 x 10^’ 

-8.38 x 10+° 

4.36 x lO- 2 

3.78 x 10+ 1 

-5.31 x 10-' 

5 

1.46 x 10 +1 

1.73 x IQ-' 

4.01 x 10+° 

5.24 x IQ- 2 

1.85 x 10 +1 

1.06 x 10-' 


TABLE 12: Jet TF parameters a, and bj for Run Ilbl and 0.4 < |rj| < 0.8. 


Light jets fojets b^ jets 


i 

Clf 

bi 

Clj 

bi 

a/ 

bi 

i 

-1.65 x 10+° 

9.52 x 10- 3 

2.67 x 10+° 

-2.07 x 10-' 

5.00 x 10 +u 

—1.15 x 10-' 

2 

3.79 x 10+° 

7.95 x 10-2 

2.79 x 10+° 

1.82 x 10-' 

1.33 x 10+° 

1.46 x 10-' 

3 

0 

2.71 x 10- 4 

0 

3.65 x 10-2 

0 

4.99 x 10- 4 

4 

1.61 x 10 +1 

-2.06 x 10-' 

-7.79 x 10+° 

3.40 x 10-2 

4.00 x 10 +1 

-8.00 x 10-' 

5 

1.47 x 10 +1 

1.72 x 10-' 

3.42 x 10+° 

6.41 x IQ- 2 

0 

4.00 x 10-' 


TABLE 13: Jet TF parameters a, and bi for Run Ilbl and 0.8 < |tj| < 1.6. 


i 

Light jets 

Clj 

bi 

b jets 

di 

bi 

bn jets 
ai 

bi 

i 

6.21 x 10+° : 

1.65 x 10-' 

7.16 x 10 +u 

-2.28 x 10-' 

8.00 x 10+° 

-1.33 x 10-' 

2 

8.73 x 10-' 

1.57 x 10-' 

4.31 x 10+° 

1.78 x 10-' 

6.04 x 10+° 

1.08 x 10-' 

3 

0 

1.77 x to -2 

0 

2.11 x 10-2 

0 

1.00 x 10- 3 

4 

2.57 x 10+° 

2.67 x 10~2 

-9.69 x 10+° 

9.76 x 10- 3 

6.00 x 10+ 1 

-7.50 x 10-' 

5 

1.01 x 10 +1 

6.36 x IQ -2 

3.41 x 10+° 

9.95 x IQ- 2 

1.00 x 10 +1 

2.00 x 10-' 


TABLE 14: Jet TF parameters a; and bi for Run Ilbl and 1.6 < |t]| < 2.5. 


Light jets b jets b ^ jets 


i 

A; 

bi 

a; 

bi 


bi 

i 

5.79 x 10+° 

-1.07 x 10-' 

1.69 x 10+° 

-2.27 x 10-' 

2.80 x 10 +1 

-3.99 x 10-' 

2 

2.98 x 10+° 

1.40 x 10-' 

3.09 x 10+° 

1.27 x 10-' 

7.99 x 10+° 

7.81 x 10-2 

3 

0 

2.61 x 10- 3 

0 

5.69 x 10- 3 

0 

7.00 x 10- 3 

4 

9.11 x 10+° 

4.16 x 10-2 

-1.08 x 10+' 

2.90 x 10-2 

9.46 x 10+° 

-9.70 x 10-2 

5 

1.72 x 10 +1 

5.85 x IQ- 2 

1.33 x 10+' 

8.71 x IQ- 2 

8.78 x 10+° 

1.13 x 10-' 




TABLE 15: 

Jet TF parameters a 

, and bi for Run IIb2 and | r\ \ < 0.4. 



Light jets 



b jets 

bn jets 


i 

Clj 

bi 

Clj 


bj dj 

bi 

1 

-2.19 x 10+° 

2.49 x 10-2 

1.22 x 10 

+0 

-1.92 x 10-' 1.60 x 10+° 

-4.93 x 10-2 

2 

4.23 x 10+° 

6.91 x 10“2 

2.55 x 10 

+0 

1.87x10-' 3.66x10+° 

8.15 x 10-2 

3 

0 

2.63 x 10- 4 

0 


4.37x10-2 0 

7.81 x 10- 3 

4 

1.41 x 10 +1 

1.81 x 10-' 

-8.57 x 10 

+0 

3.83 x 10-2 7 .48 x 10+° 

-1.91 x 10-' 

5 

1.48 x 10 +1 

1.80 x 10-' 

3.85 x 10 

+0 

5.67x 10-2 4.56x 10+° 

1.59 x 10-' 



39 


TABLE 16: Jet TF parameters a; and bj for Run IIb2 and 0.4 < |rj| < 0.8. 


i 

Light jets 

a, bj 

b jets 

cij bj 

bfi jets 

aj bj 

1 

-1.68 x 10 +u 

1.37 x 10~ 2 

2.19 x 10+ u 

-2.02 x 10- 1 

6.00 x 10 +u 

-1.46 x 10~' 

2 

4.12 x 10+° 

7.83 x 10~ 2 

2.53 x 10+° 

1.86 x 10- 1 

1.38 x 10+° 

1.60 x 10~’ 

3 

0 

2.80 x 10~ 4 

0 

3.43 x 10~ 2 

0 

1.50 x 10~ 4 

4 

1.68 x 10 +1 

-2.16 x 10^’ 

-8.54 x 10+° 

3.46 x 10~ 2 

3.75 x 10+ 1 

-3.59 x 10-' 

5 

1.56 x 10 +1 

1.66 x 10~' 

3.59 x 10+° 

6.47 x 10~ 2 

1.59 x 10 +1 

1.87 x 10~' 


TABLE 17: Jet TF parameters a, and bj for Run IIb2 and 0.8 < | rj j < 1.6. 


i 

Light jets 

aj bj 

b jets 

ctj bj 

bn jets 

cij bj 

1 

2.20 x 10 +u 

-2.52 x 10^ 2 

6.75 x 10+ u 

-2.27 x 10- 1 

8.00 x 10 +u 

— 1.39 x lO^ 1 

2 

5.77 x 10+° 

1.07 x 10~’ 

4.14 x 10+° 

1.78 x 10- 1 

5.54 x 10+° 

1.28 x 10“’ 

3 

0 

3.57 x 10~ 4 

0 

1.93 x 10~ 2 

0 

1.43 x 10~ 4 

4 

2.16 x 10 +1 

-1.23 x 10~' 

-1.09 x 10+ 1 

1.23 x 10~ 2 

—4.00 x 10 +1 

-2.22 x 10 _1 

5 

1.91 x 10 +1 

1.19 x 10~' 

3.25 x 10+° 

1.03 x lO^ 1 

4.10 x 10 +1 

2.87 x 10 _1 


TABLE 18: Jet TF parameters a, and bj for Run IIb2 and 1.6 < |tj| < 2.5. 


i 

Light jets 

at bj 

b jets 

cij bj 

bn jets 

cij bj 

1 

7.33 x 10 +IJ 

-1.41 x 10 _l 

2.55 x 10+ u 

-2.48 x 10- 1 

9.25 x 10 +u 

-1.81 x lO^ 1 

2 

3.16 x 10+° 

1.35 x 10 _1 

4.56 x 10+° 

1.13 x 10 _1 

3.35 x 10+° 

1.53 x 10 _1 

3 

0 

3.51 x 10~ 3 

0 

5.90 x 10~ 3 

0 

1.68 x 10~ 4 

4 

5.54 x 10+° 

4.50 x 10~ 2 

-1.21 x 10+' 

2.30 x 10~ 2 

2.40 x 10 +1 

2.36 x 10~' 

5 

1.60 x 10 +1 

6.68 x 10~ 2 

1.32 x 10+' 

8.86 x 10~ 2 

3.55 x 10+ 1 

—1.18 x 10 _1 


TABLE 19: Jet TF parameters a,- and b t for Run IIb3 and |tj| < 0.4. 


i 

Light jets 

cij bj 

b jets 

ctj bj 

bfj jets 

at bj 

1 

-1.87 x 10 +u 

1.94 x 10~ 2 

4.64 x 10 +u 

-2.02 x lO^ 1 

6.72 x 10 +u 

-1.39 x 10 _1 

2 

4.11 x 10+° 

7.20 x 10~ 2 

2.78 x 10+° 

2.32 x 10 _1 

1.56 x 10 +0 

1.53 x 10~' 

3 

0 

3.39 x 10~ 4 

0 

7.93 x 10~ 2 

0 

2.60 x 10~ 4 

4 

1.70 x 10 +1 

-1.47 x 10~' 

-6.31 x 10+° 

1.41 x 10~ 2 

3.27 x 10 +1 

-3.09 x 10~' 

5 

1.63 x 10 +1 

1.60 x 10~’ 

3.56 x 10+° 

6.45 x 10~ 2 

1.32 x 10 +1 

1.98 x 10~' 


TABLE 20: Jet TF parameters a, and bj for Run IIb3 and 0.4 < 1 77 1 < 0.8. 


i 

Light jets 

cij bj 

b jets 

at bj 

bn jets 

cij bj 

1 

— 1.70 x 10+° 

1.67 x 10~ 2 

4.22 x 10 +o 

-2.03 x lO^ 1 

6.70 x 10+° 

-1.44 x 10 _1 

2 

4.17 x 10+° 

7.83 x 10~ 2 

2.91 x 10+° 

2.08 x 10 _1 

1.64 x 10 +0 

1.56 x 10~' 

3 

0 

3.12 x 10~ 4 

0 

4.86 x 10~ 2 

0 

1.98 x 10~ 4 

4 

1.89 x 10 +1 

— 1.77 x 10~’ 

-7.71 x 10° 

2.99 x 10~ 2 

3.68 x 10 +1 

-2.84 x 10~' 

5 

1.77 x 10 +1 

1.55 x 10~' 

3.43 x 10+° 

6.77 x 10~ 2 

1.66 x 10 +1 

1.76 x 10~' 
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TABLE 21: Jet TF parameters a; and bj for Run IIb3 and 0.8 < |jj| < 1.6. 


i 

Light jets 

fl; bj 

b jets 

at bi 

bn jets 

a, b. 

1 

1.56 x 10 +u 

-1.84 x 10~ 2 

1.02 x 10+' 

-2.61 x 10- 1 

8.00 x 10 +u 

-1.33 x 10 1 

2 

5.96 x 10+° 

1.08 x 10“’ 

4.92 x 10+° 

1.91 x 10- 1 

6.04 x 10+° 

1.08 x 10“’ 

3 

0 

3.87 x 10~ 4 

0 

3.34 x 10~ 2 

0 

1.00 x 10~ 3 

4 

2.12 x 10 +1 

-9.65 x 10~ 2 

-9.51 x 10+° 

-5.50 x 10- 3 

6.00 x 10+ 1 

-7.50 x 10“’ 

5 

2.06 x 10 +1 

1.08 x 10 _1 

3.00 x 10+° 

1.16 x 10 _1 

1.00 x 10 +1 

2.00 x 10 _1 


TABLE 22: Jet TF parameters a; and bj for Run IIb3 and 1.6 < 1 77 1 < 2.5. 


i 

Light jets 

at bj 

b jets 
a-, 

bi 

bfj jets 

«i bj 

1 

6.59 x 10 +u 

— 1.06 x 10“ 1 

-3.29 x 10 1 

— 1.80 x lO^ 1 

2.35 x 10 +I 

-3.25 x lO^ 1 

2 

2.77 x 10+° 

1.37 x 10 _1 

9.64 x 10~ ! 

1.47 x lO^ 1 

6.73 x 10+° 

1.00 x 10 _1 

3 

0 

2.01 x 10~ 3 

0 

3.62 x 10~ 3 

0 

3.65 x 10~ 3 

4 

1.26 x 10 +1 

2.99 x 10~ 2 

—6.32 x 10 +0 

1.59 x 10~ 2 

9.75 x 10+° 

-6.83 x 10~ 2 

5 

1.78 x 10 +1 

5.45 x 10~ 2 

1.47 x 10 +1 

7.99 x 10~ 2 

1.29 x 10 +1 

9.45 x 10~ 2 



